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1.  INTRODUCTION  AND  SUMMARY 


Closed  cycle  operation  of  CO2  electric  discharge  lasers  has  been  shown 
to  be  substantially  influenced  by  the  presence  of  molecular  species  with  high 
electron  attachment  cross  sections,  negative  ions  which  increase  the  plasma 
heating,  and  species  which  efficiently  deactivate  the  upper  laser  level  (Refs. 

1  and  2).  Depending  on  the  origin  of  these  electrophilic  species  we  can  dis¬ 
tinguish  basically  three  kinds.  First,  impurities  which  are  certain  trace 
species  present  in  the  initial  laser  gas  mixture  before  the  discharge.  Second, 
outgassing  and  material  degradation  of  the  device  hardware  produce  what  we 
shall  call  contaminants.  Finally,  plasma  by-products  are  produced  from  the 
ideal  laser  gas  mixture  during  and  after  the  discharge  due  to  plasma  chemistry. 
The  present  study  is  mainly  concerned  with  plasma  by-products  and,  to  some 
degree,  with  initial  impurities  and  their  influence  on  laser  performance. 

The  plasma  chemistry  important  in  the  formation  of  these  by-products 
has  been  studied  in  greatest  detail  for  He/N^/CO^  mixtures  loaded  by  a  dc 
discharge.  The  Army  and  the  Air  Force  are  currently  pursuing  methods  for 
eliminating  helium  mixtures  in  repetitively  pulsed  laser  systems.  Thus, 
particularly  for  closed  cycle  operation,  there  exists  interest  in  determining 
dominant  processes  detrimental  to  sustained  repetitive  high  energy  density 
pulsed  loading  of  generic  I^/^/CO^  laser  mixtures  which  might  also  con¬ 
tain  CO,  H2O  or  other  species  characteristic  of  hydrocarbon  fuel  combustion 
products. 

The  purpose  of  the  effort  described  here  was  to  identify  chemical  re¬ 
actions  leading  to  the  generation  of  plasma  by-products  in  several  laser  gas 
mixtures,  to  determine  the  relative  importance  of  the  various  plasma  by-products 
in  degrading  assessed  laser  performance,  and  to  furnish  predictions  of  all  im¬ 
portant  molecular  species  production  and  decay  rates  during  and  following 


the  electrical  pulse.  In  doing  this,  maximum  use  was  to  be  made  of  existing 
kinetic  data  and  of  the  existing  continuous  wave  electric  laser  code  previously 
developed  by  the  authors  (Ref.  3). 

Furthermore,  the  steady  state  buildup  of  byproduct  concentrations  under 
repetitively  pulsed  operations  was  to  be  examined,  and  key  experimental  efforts 
to  verify  predicted  results  were  to  be  recommended. 

The  essential  features  of  the  modeling  approach  used  in  this  effort  are 
described  in  Section  2.  The  discussion  includes  the  basic  equations  as  well 
as  some  details  of  the  plasma  kinetics  model  employed.  Results  obtained  for 
several  laser  gas  mixtures  are  described  and  compared  with  experimental 
data  in  Section  3.  The  good  agreement  achieved  lends  credibility  to  the  theo¬ 
retical  model  and  suggests  that  application  of  the  analysis  to  other  laser  gas 
mixtures  may  be  made  with  confidence. 

Section  4  presents  a  detailed  analysis  of  test  cases  for  two  closed  cycle 
EDL.  devices  currently  under  development.  The  study  includes  the  effects  on 
performance  of  variations  in  the  electric  field,  the  effective  electron  beam 
current,  gas  mixture  composition,  cycle  simulation  and  gas  contamination  by 
water  vapor. 

Section  5  presents  the  conclusions  derived  from  the  work  accomplished. 
Also  presented  are  recommendations  for  key  experimental  efforts  that  should 
be  undertaken  to  verify  predicted  results  as  well  as  for  future  theoretical 
efforts  that  should  be  performed  to  improve  the  computer  model. 

Briefly  summarizing,  a  theoretical  model  is  presented  which  can  be 
used  to  analyze  plasma  chemistry  effects  on  the  performance  of  pulsed,  elec¬ 
tron  beam  sustained  electric  discharge  lasers.  The  one-dimensional  model 
contains  descriptions  of  the  configuration  components  of  primary  importance 
which  are  the  electron  beam,  the  discharge,  the  plasma  and  its  composition, 
laser  radiation,  and  their  mutual  interaction.  All  important  processes  are 
modeled  in  terms  of  experimentally  controllable  parameters  which  include 


the  post-foil  electron  beam  current  density,  the  primary  electron  energy,  the 
applied  electric  field  potential,  discharge  electrode  spacing,  initial  laser  gas 
composition,  temperature  and  pressure,  both  electron  beam  and  electric  dis¬ 
charge  pulse  widths  and  relative  sequencing,  as  well  as  cavity  geometric  data 
such  as  mirror  reflectivities,  transmissivities  and  gain  length. 

Finally,  a  complete  list  of  rate  mechanisms  and  coefficients  as  used  in 
this  work  is  presented  in  an  Appendix  to  this  report. 
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2.  MODELING  APPROACH 

2.1  CAVITY  DISCHARGE  MODEL 

This  section  presents  a  brief  description  of  the  essential  features  which 
a  theoretical  model  must  possess  to  adequately  describe  plasma  chemistry 
processes  in  pulsed  electric  discharge  lasers.  This  model  incorporates  de¬ 
scriptions  of  the  electron  beam,  the  discharge,  the  plasma  composition  and 
laser  radiation,  and  their  mutual  interaction.  The  effect  of  both  single  and 
repetitive  pulsing  can  be  modeled.  To  simulate  selective  removal  of  detri¬ 
mental  species  (scrubbing)  the  plasma  composition  can  be  modified  between 
successive  pulses  when  recirculation  and  repetitive  pulsing  are  to  be  con¬ 
sidered.  Modeling  as  many  of  the  basic  features  of  a  device  as  possible  serves 
to  minimize  simplifying  assumptions  which  always  have  to  be  made.  One  such 
simplification  made  in  this  work  is  the  assumption  of  spatially  uniform  con¬ 
ditions  in  the  cavity  resulting  in  a  one -dimensional  model.  On  the  other  hand, 
this  analysis  fully  counles  the  energy  equation  into  the  system  of  species  rate 
equations,  eliminating  the  often  made  assumption  that  variations  in  the  gas 
temperature,  which  is  of  prime  importance  to  laser  operation,  can  be  ignored. 

Considering  uniformly  excited  laser  cavity  volumes  of  the  order  of  a  few 
liters,  and  electric  discharge  pulses  of  approximately  microsecond  duration  it 
can  easily  be  shown  that  a  constant  density,  time-dependent  treatment  of  the 
plasma  chemistry  should  adequately  describe  the  important  processes  occur¬ 
ring  in  a  pulsed  discharge.  For  repetitively  pulsed  discharges  and  closed  loop 
operation  of  the  laser,  the  flow  velocity  around  the  cycle  is  assumed  to  be  suf¬ 
ficiently  small  that  pressure  changes  due  to  varying  flow  velocity  can  be  neg¬ 
lected.  Plasma  composition  changes  around  the  recirculation  loop  can  then  be 
treated  in  the  same  time-dependent  manner  as  in  the  cavity,  provided  that  post¬ 
discharge  expansion  of  the  plasma  as  well  as  cooling  by  heat  exchangers  is 
appropriately  accounted  for. 
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The  laser  gas  mixture  consists  of  neutral  and  charged  particles,  vibra- 
tionally  and  electronically  excited  particles,  and  electrons.  Therefore  a 
multitude  of  species  must  be  tracked  during  the  calculations.  It  is  assumed 
that  different  quantum  states  of  various  particles  can  be  treated  as  separate 
species.  For  simplicity,  it  is  also  assumed  that  a  single  vibrational-rotational 
radiative  transition  sufficiently  accounts  for  laser  output.  While  not  generally 
true,  this  is  a  reasonable  assumption  if  lasing  of  the  CO^  molecule  is  considered. 

The  model  outlined  in  this  section  adequately  represents  phenomena 
occurring  in  the  positive  column.  While  a  detailed  description  of  cathode 
sheath  phenomena  is  outside  the  scope  of  this  work,  a  procedure  will  be 
outlined  at  the  end  of  this  section  which  permits  us  to  estimate  the  effect  of 
heating  in  the  cathode  sheath  on  plasma  chemistry  and  the  generation  of  plasma 
byproducts. 

Under  the  above  stated  assumptions  the  basic  governing  equations  can 
be  written  as  follows: 


Continuity: 


_ 

dt  " 


0 


Species: 


dF. 

f>  dT  =  *i  +  <6il  -6i2^ 


(1) 

(2) 


Here  p  denotes  the  density,  F^  the  mol/mass  ratio  of  the  i*^  species,  and  w. 
its  overall  net  rate  of  production  due  to  any  reaction  mechanism  except  for 
radiation.  For  an  arbitrary  reaction  of  the  form 


N  N 

^  y;  Ai^r  £  v{'  A. 


i  =  1 


i  =  l 


(3) 
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the  net  rate  of  production  for  any  participating  species  for  which  p."  -  pi  /  0 
can  be  written  as 

W.  =  (pi'  -  v!)  jkf  n  (pF.)  -  kb_/7  (pF.)v  j  (4) 

where  k^  and  k^  denote  the  forward  and  backward  rate  constant,  respectively. 
By  the  law  of  mass  action,  this  formulation  presumes  thermal  equilibrium  of 
the  species  at  a  temperature,  T  (or  electron  energy,  u).  The  rate  constants 
are  then  functions  of  T  (or  u). 


Laser  radiation,  its  intensity  denoted  by  I,  depletes  the  upper  laser  level 
and  replenishes  the  lower  laser  level,  as  indicated  by  the  Kronecker  symbols 
in  Eq.  (2).  Gain  and  photon  energy  are  denoted  by  a  and  e  ,  respectively. 

Energy: 

i(pe  )--jDE-,£i  (5) 

which  simply  states  that  any  change  of  internal  energy  of  the  laser  gas  is  due 
to  discharge  energy  pumped  into  the  plasma,  or  due  to  radiation  leaving  the 
gas  in  the  cavity.  Here  denotes  the  drift  field  current  density  and  E  the 
electric  field  strength.  The  quantity  K  represents  a  correction  factor  and  is 
defined  as 

*  g  i.-_g*pL-gy  (6) 

aL 


where  a  is  the  cavity  threshold  gain.  The  set  of  governing  equations  is  com¬ 
pleted  by  adding  an  equation  of  state.  Species  thermodynamic  data  including 
specific  heats,  entropy  and  enthalpy  are  used  in  tabular  form  for  all  species 
considered. 


Assuming  a  Fabry-Perot  cavity,  the  threshold  gain  can  be  expressed  in 
terms  of  mirror  reflectivities  and  optical  path  length  through  the  gain  medium 
as 


a 


fn(rlr2 

2L 


(7) 


6 


Note  that  for  aL  «  1,  which  is  most  often  the  case,  k  reduces  to  unity.  Lasing 
commences  when  the  medium  gain  exceeds  the  cavity  gain,  and  ceases  when  the 
pumping  reactions  fail  to  maintain  the  medium  gain  at  the  cavity  gain  level. 

The  medium  gain  can  be  expressed  as 

a  =  pe(S21r2-S12r1)V  (8) 

where  end  F ^  are  the  mol/mass  ratios  of  the  upper  and  the  lower  laser 
level,  respectively.  All  spectroscopic  data  enter  the  equation  through 
and  S  ^2  which  are  functions  of  temperature  only.  V  denotes  the  Voigt  function 
which  basically  represents  a  generalized  line  shape  function, 

V(rj)  =  exp(n2)  [1  -  erf(rj)]  (9) 


with  q  being  the  ratio  of  the  Lorentz  line  width  to  the  Doppler  line  width,  i.e.. 


n 


(ini) 


1/2  A"L 


D 


(10) 


The  set  of  basic  equations  is  completed  by  the  equation  of  state, 

NS 

P  =  PS  T  2  Fi  (11) 

i=  1 


which  determines  the  pressure  as  a  function  of  density,  temperature  and  gas 
composition. 


Energizing  the  laser  gas  in  an  electric  discharge  requires  the  presence 
of  free  electrons  which  are,  at  least  initially,  here  assumed  to  be  generated 
by  electron-beam  ionization.  Assuming  that  both  the  effective  electron-beam 
current  density,  jj-g»  and  the  effective  electron  beam  energy,  Egg*  are  speci¬ 
fied,  the  number  density  of  primary  electrons  in  the  laser  gas  can  be  expressed 
as 


JEB 

EB  =  %  VEB 


(12) 
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where  is  the  electron  charge,  and  VEB»  the  relativistic  electron  beam 
velocity,  is  given  by 


V 


EB 


(13) 


Here  c  denotes  the  speed  of  light,  and  ED-,„  the  electron  rest 
equivalence  (mec  =  511  keV).  Knowing  the  primary  electron 
rate  of  ionization  and  production  of  secondary  electrons  due  to 
beam  only  can  be  determined  according  to  the  reaction 


mass  energy 
density,  the 
the  electron 


Ai  ♦  "EB 


£,  i  A+  . 

Ai  +  ne  +  nEB 


(14) 


where  the  primary  ionization  rate  constant,  k£  ^  can  be  evaluated  as  a  func¬ 
tion  of  the  ion  mass,  the  beam  ionizing  power  and  the  beam  electron  velocity 
(Section  2.2). 

The  drift  field  current  density,  and  therefore  the  energy  which  can  be 
pumped  into  the  gas,  is  proportional  to  the  number  of  free  electrons  and  their 
drift  velocity  such  that,  to  good  approximation, 

jD  =  ne%VD  <15) 

The  drift  velocity,  Vp,  is  a  function  of  the  gas  mixture  composition  and  the 
electric  field/density  ratio.  It  is  obtained  from  an  evaluation  of  the  electron 
energy  distribution  function  (Boltzmann  code,  see  Section  2.2)  and  can  be 
expressed  as 


(16) 
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The  volumetric  rate  of  energy  addition  in  the  discharge  can  then  be  expressed 
as 


n  q 
e  e 


r\n 

C1 +  c2(l 

!)j 

(17) 


where  C^,  C ^  and  n  are  mixture  dependent  constants. 


Equation  (17)  is  not  only  used  to  evaluate  the  energy  added  in  the  dis¬ 
charge  but  also  to  determine  the  secondary  electron  energy  which  is  essential 
in  evaluating  electron  impact  excitation  rates.  Postulating  a  basic  electron 
impact  excitation  mechanism  as 

k,  .  ^ 

Ai  *  e— '-L  Al  +  e  (18) 


the  power  consumed  by  this  mechanism  can  be  expressed  as 

AP  =  L  nenikf,i  ^i  (19) 

i 

where  the  summation  indicates  that  more  than  one  mechanism  may  be  involved. 
In  fact,  as  will  be  seen  later  (see  Section  2.2)  the  summation  extends  over  exci¬ 
tation  mechanisms  including  vibrational  excitation,  electronic  excitation,  and 
Ionization.  The  k,  .  which  appear  here  are  those  which  normally  result  from 
a  full  Boltzmann  analysis  of  electron  kinetics  in  the  primary  laser  gas  but 
exclusive  of  electron  kinetics  in  the  plasma  discharge  byproducts. 


Applying  conservation  of  energy  by  equating  the  right-hand  sides  of 
Eqs.  (17)  and  (19)  then  leads  to  the  important  relation 

%E[CltC2(l)n]=  <20> 

Since  the  impact  excitation  rate  constants,  k,  .  can  be  expressed  as  functions 
of  the  electron  energy  which  itself  is  a  function  of  E/N,  Eq.  (20)  constitutes 
an  implicit  relation  for  the  electron  energy,  and  simultaneously  states  its 
functional  dependence  on  E/N  in  terms  of  input  functions.  It  is  through  Eq.  (20) 
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i 


**  - 


that  coupling  of  the  Boltzmann  code  and  the  discharge  calculations  is  put  into 
effect.  If  all  electron  excitation  mechanisms  considered  in  the  Boltzmann  code 
are  included  on  the  right-hand  aide  of  Eq.  (20),  the  average  electron  energy 
iteratively  computed  from  Eq.  (20)  will  accurately  compare  with  the  average 
electron  energy  computed  by  the  Boltzmann  code  for  the  given  E/N. 

The  reason  for  using  the  Boltzmann  code  in  terms  of  Eq.  (20)  is  both 
historical  and  practical.  In  earlier  work  we  did  not  have  a  Boltzmann  code. 

We  therefore  had  to  fit  available  results  which  then  led  to  Eq.  (20).  The 
practical  reason  for  retaining  this  approach  is  the  fact  that  the  Boltzmann 
code  requires  more  than  twice  the  core  storage  used  by  the  remainder  of 
the  analysis  rendering  inclusion  of  the  Boltzmann  code  as  a  subroutine  both 
costly  and  impractical. 

Equations  (1)  through  (20),  in  conjunction  with  a  detailed  set  of  reactions 
and  rate  coefficients,  to  be  discussed  in  Section  2.2,  provide  the  basic  frame¬ 
work  to  evaluate  plasma  conditions  in  the  positive  column  of  the  discharge 
as  a  function  of  time,  and  if  desired,  the  radiative  output  of  the  laser  cavity. 

In  the  remainder  of  this  subsection  an  approach  will  be  outlined  to  assess 
plasma  chemistry  effects  in  the  electrode  sheaths.  Because  the  heating  in  the 
anode  sheath  is  not  intense  (Ref.  2)  we  can  confine  our  considerations  to  the 
cathode  fall  region. 

In  the  thin  cathode  sheath  region,  the  local  heating  rate  is  very  high. 

This  produces  a  relatively  high  local  pressure  initially,  and  interferograms 
indicate  that  changes  of  density  are  small  initially  (Ref.  4).  Owing  to  the 
high  temperatures  which  accompany  heating  near  the  cathode  it  can  be  ex¬ 
pected  that  heat  transfer  from  the  gas  to  the  cathode  cannot  be  ignored.  On 
the  other  hand,  as  pointed  out  by  Long  (Ref.  5)  convective  cooling  of  the 
cathode  sheath  can  be  ignored  if  the  flow  velocity  in  the  positive  column  is 
small  enough.  Ignoring  convective  cooling  therefore  is  consistent  with  our 
model  presented  here.  In  order  to  compute  the  gas  composition  in  the  sheath, 
knowledge  of  the  temperature  in  that  region  is  most  important.  According  to 


•at 


mm* 


Ref.  5,  the  maximum  temperature  rise  in  the  cathode  sheath  as  a  function  of 
bulk  energy  input  and  heat  transfer  to  the  cathode  can  be  estimated  from 


To  = 


1  c 

5  ir~ 


(21) 


where  T  is  the  maximum  temperature  in  die  sheath,  T  is  the  cathode  wall 
c  o 

temperature,  kQ  is  the  thermal  conductivity  of  the  gas,  and  dc  is  the  cathode 
sheath  thickness.  Q  represents  the  bulk  energy  input  into  the  cathode  sheath. 
Assuming  current  continuity,  Q  is  given  by 

Q  =  JD  Ec  (22) 

where  is  the  discharge  current  density,  and  E  represents  the  average 
electric  field  in  the  cathode  sheath. 


To  compute  the  temperature  increase  according  to  Eq.  (21)  the  cathode 
sheath  thickness  and  the  cathode  fall  potential  must  be  known.  Values  for 
these  quantities  are,  unfortunately  extremely  scarce.  Since  the  dominant 
constituent  of  laser  gases  of  interest  here  is  N^,  and  since  we  can  only  ex¬ 
pect  to  obtain  some  more  or  less  qualitative  information  using  this  simplified 
approach  anyhow,  results  presented  in  the  form  of  Fig.  1  for  will  be  used 
to  obtain  values  for  the  cathode  fall  potential  Vc>  In  Fig.  1,  pQ  denotes  the 
reduced  pressure  at  the  cathode  surface,  i.e., 


=  P 


273(H) 


(23) 


The  procedure  described  above  yields  a  maximum  temperature  in  the  cathode 
sheath  which  is  in  excellent  agreement  with  the  example  given  in  Ref.  5. 


Having  thus  determined  a  representative  cathode  sheath  temperature, 
the  initial  pressure  is  given  by  the  equation  of  state,  assuming  constant  den¬ 
sity.  It  is  reasonable  to  assume  that  the  relatively  high  temperature  and 
pressure  will  promote  quick  establishment  of  thermochemical  equilibrium. 
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Fig.  1  -  Voltage -Current  Characteristic  and  Cathode  Fall  Distance  in  N2  (Ref.  5) 


Therefore,  under  these  conditions,  the  resulting  gas  composition  is  easily 
computed  using  the  CEC  program  (Ref.  6). 


The  overall  effect  of  possible  contaminants  generated  in  the  thin  cathode 

sheath  on  the  plasma  chemistry  processes  occurring  in  the  positive  column, 

which  forms  the  bulk  of  the  discharge,  must  take  into  consideration  the  different 

conditions  involved.  Resulting  concentrations,  averaged  over  the  sheath  and  the 

positive  column,  can  be  obtained  by  performing  a  non-reactive  mixing  calcula- 

th 

tion,  which  defines  the  average  mole  fraction  of  the  i  species  as 


X,  = 


n.  +  n. 
■l?-P _ ii£_ 


1  y  n.  tVn. 


(24) 


th 

where  n.  and  n.  are  the  number  densities  of  the  iw  species  in  the  positive 
i.p  l,  c  r  r 

column  and  in  the  cathode  sheath,  respectively.  Expressing  the  number  den¬ 
sities  in  terms  of  pressure,  temperature,  volume  and  mole  fraction,  Eq.  (24) 
can  also  be  written  in  the  more  convenient  form 


Xi..»tXi.c 


R 


1  +  R 


(25) 


where 


(26) 


Because  the  volume  of  the  cathode  sheath  is  likely  to  be  orders  of  magnitude 
smaller  than  that  of  the  positive  column,  and  because  T£  »  Tp  due  to  heating 
of  cathode  gases,  we  can  reasonably  expect  that  R  «  1,  and  then  Eq.  (*5)  shows 
that  concentrations  of  plasma  byproducts  formed  in  the  cathode  fall  region  will 
have  to  be  extremely  large  to  have  a  measurable  effect  on  the  overall  gas  com¬ 
position. 
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2.2  PLASMA  KINETICS 
2.2.1  Primary  Ionization 

In  order  to  obtain  high  efficiency  in  electric  discharge  lasers  it  has  been 
found  that  it  is  advantageous  to  separate  excitation  into  primary  and  secondary 
processes,  each  tailored  to  its  task.  Primary  ionization  can  be  realized  by 
various  means,  most  currently  using  an  electron  beam. 

To  predict  the  yield  of  secondary  electrons  due  to  electron  beam  ioniza¬ 
tion  according  to  Eq.  (14)  we  need  to  know  k,  . ,  the  primary  ionization  rate  con- 
stant.  The  rate  at  which  secondary  electrons  are  produced  may  be  written  as 


where  p.  is  the  mass  density  of  the  ith  species,  E^  is  its  effective  ionization 
energy,  and  (dE/dm).  is  the  primary  electron  energy  loss  rate  due  to  interaction 
with  the  ith  species  of  the  gas. 

From  Eq.(14)  we  obtain  similarly 


dn 

e 

"3T  = 


(28) 


Combining  the  above  two  expressions  ana  using  Eq.(12)  it  is  found  that 


c  /Wi\  VEB  .  dE  . 
'£*L  ~  NA  \Ei/  C  dm  i 


(29) 


Figure  2  shows  V__  and  (dE/dm).  for  various  species  of  interest  as  a  func- 
tion  of  effective  electron  beam  energy.  Noticing  that  the  product  of  Vgg  and 
(dE/dm).  is  a  function  of  Egg  as  shown  in  Fig.  3,  the  primary  ionization  rate 
constant  ^  can  be  expressed  as  a  function  of  Egg  only.  Using  the  curves 
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shown  in  Fig.  3  it  is  found  that  an  EEE  dependence  reproduces  those  curves 
to  within  5%  for  40  <  E__  (keV)  <  300. 

2.2.2  Secondary  Excitation 

In  order  to  analyze  in  quantitative  detail  the  electron,  atomic  and/or 
molecular  physics  of  electrically  excited  laser  discharges,  it  is  necessary 
to  accurately  determine  all  of  the  laser  level  excitation  rates  and  transport 
properties  which  determine  the  gross  properties  of  these  discharges  (Ref.  8). 
To  obtain  these  parameters,  knowledge  of  the  electron  distribution  function  is 
required.  Early  work  assumed  that  the  distribution  was  Maxwellian,  but  more 
recent  work  has  shown  this  assumption  to  be  incorrect  (Ref.  9).  Consequently, 
it  is  necessary  to  numerically  solve  the  Boltzmann  transport  equation  for  the 
particular  value  of  E/N  and  gas  mixture  in  order  to  obtain  accurate  values  for 
these  excitation  rates  and  transport  properties. 

The  Boltzmann  code  used  in  the  present  work  was  obtained  from  the 
Air  Force  Weapons  Laboratory  (Ref.  10).  The  code  uses  the  method  of  finite 
differences  to  solve  the  Boltzmann  equation  for  a  spatially  uniform  gas  in 
the  presence  of  a  steady  electric  field,  as  described  by  Nighan  (Ref.  9). 

Elastic  heating,  vibrational  excitation,  ionization,  electronic  excitation  and 
disociative  attachment  processes  are  included,  however,  superelastic  and 
electron-electron  interactions  are  not.  The  electron  impact  cross  sections 
employed  here  are  the  same  as  computed  by  Phelps  and  coworkers  (Ref.  11). 

The  output  data  consist  of  the  fractional  discharge  energy  transfer  to 
the  various  excitation  modes  and  the  corresponding  rate  constants  as  well 
as  several  macroscopic  parameters  such  as  electron  drift  velocity,  char¬ 
acteristic  energy,  diffusion  coefficient,  and  effective  average  electron  energy. 
The  latter  is  defined  such  that,  for  a  Maxwellian  energy  distribution,  it  corre¬ 
sponds  to  the  electron  temperature,  i.e., 

00 

ug  =  J'  u3/2  f  (u)  du  (30) 
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where  u  denotes  electron  energy  and  f(u)  is  the  distribution  function.  An 
energy  conservation  check  is  performed  at  the  end  of  each  computation  to 
indicate  the  validity  of  the  results. 

This  code  has  been  designated  primarily  for  use  with  CO^  mixtures. 
Angularly  averaged  cross-section  data  are  provided  for  CO^,  N^,  CO, 

H  ,,  Oz  and  He.  The  vibrational  levels  have  been  separately  considered  but 
several  electronic  and  ionization  levels  for  each  gas  have  been  lumped  to¬ 
gether.  Addit  onal  modifications  have  been  incorporated  in  order  to 
obtain  consistency  between  the  Boltzmann  code  and  the  general  kinetic 
rate  model  used  in  the  present  work.  Specifically,  to  keep  the  number 
of  reaction  mechanisms  a,nd  species  manageable,  the  consideration  of  CO, 

Ct 

vibrational  levels  is  limited  to  00^(000),  00^(100),  CO^OnO  with  n  =  1,2,3 
and  00^(001).  For  all  diatomic  molecules  only  the  first  vibrational  level  is 
considered  in  the  plasma  chemistry  aside  from  the  ground  level.  Since  the 
Boltzmann  code  considers  up  to  eight  vibrational  levels  for  each  molecule, 
simply  ignoring  the  upper  levels  would  have  led  to  an  inconsistency  in  the 
energy  balance.  Instead,  we  added  an  additional  calculation  which  computes 
from  the  rate  constants  for  the  upper  vibrational  levels  an  energy  weighted 
average  rate  constant  consistent  with  conservation  of  energy  and  the  reaction 
mechanisms  used  in  the  plasma  chemistry  analysis.  The  number  of  levels 
over  which  the  averaging  is  performed  can  be  selected  by  the  program  user 
depending  on  how  many  of  the  higher  vibrational  levels  of  the  same  molecule 
are  considered  in  the  plasma  kinetics. 

Transport  data  consisting  of  the  drift  velocity  and  the  average  electron 
energy  for  various  laser  gas  mixtures  are  shown  in  Figs,  4  and  5,  respectively. 
For  a  given  electron  density,  the  drift  velocity  determines  the  discharge  cur¬ 
rent  density  and  thus,  when  multiplied  by  the  electric  field  strength,  determines 
the  volumetric  energy  addition  to  the  gas  in  the  discharge.  For  each  gas  mix¬ 
ture  the  drift  velocity  shown  in  Fig.  4  is  curve  fitted  according  to  Eq.  ( 1 6)  for 
use  in  computing  the  average  electron  energy  from  Eq.(20).  The  drift  velocity 
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Fig.  4  -  Electron  Drift  Velocity  as  a  Function  of  E/N 
for  Various  Gas  Mixtures 
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Fin. 5  -  Average  Electron  Energy  as  Function  of  E/N  for  Various 
Gas  Mixtures 
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provides  one  of  the  links  that  couple  the  Boltzmann  code  to  the  plasma  chem¬ 
istry  analysis,  the  other  being  the  impact  excitation  rates.  Drift  velocity  curve 
fit  coefficients  for  various  gas  mixtures  are  shown  in  Table  1.  These  curve  fits 
are  very  accurate  over  a  wide  range  of  E/N  (1  •  10"^<  E/N  (Vcm^)<  7  •  10~*^). 
The  relationship  between  average  electron  energy  and  E/N  shown  in  Fig.  5  is 
not  used  directly  in  the  analysis  but  is  useful  for  checking  purposes.  If  all 
relevant  impact  excitation  mechanisms  have  been  included  in  determining  the 
average  electron  energy  from  Eq.  (20),  the  value  obtained  will  closely  agree 
with  the  value  computed  by  the  Boltzmann  code,  and  shown  in  Fig.  5  as  a  func¬ 
tion  of  E/N  and  gas  mixture. 

Fractional  discharge  energy  transfer  into  various  excitation  levels  is 
shown  in  Figs.  6  through  9  for  a  number  of  different  gas  mixtures.  Although 
the  plasma  chemistry  analysis  makes  no  direct  use  of  this  information,  it  is 
illustrative  as  well  as  useful  in  deciding  at  which  value  of  E/N  a  given  discharge 
would  be  most  efficient.  For  this  purpose  we  have  shown  on  each  figure  that 
fraction  of  the  discharge  energy  which  is  transferred  into  the  upper  laser  level 
and  into  other  excitation  levels  which  in  turn  transfer  their  energy  to  the  upper 
laser  level.  The  location  of  the  peak  should  be  a  good  indication  of  the  optimum 
E/N  value  for  the  discharge.  It  is  interesting  to  note  that  the  nominal  oper¬ 
ating  value  of  E/N  =  3.7  •  10  ^  Vcm^  of  the  ABEL  device  (analysis  of  which  is 
discussed  in  Section  3)  is  almost  the  exact  location  of  the  maximum  of  the  re¬ 
spective  curve  in  Fig.  8.  It  is  also  observed  that  mixtures  containing  appear 
to  show  an  optimum  value  of  E/N  which  is  distinctly  higher  than  that  for  mix¬ 
tures  containing  He.  Rate  constants  for  the  various  electron  impact  excitation 
mechanisms  forming  the  other  link  between  the  Boltzmann  code  and  die  plasma 
chemistry  and  laser  calculations  are  shown  in  Figs.  10  through  20.  These  mech¬ 
anisms  include  vibrational  excitation,  electronic  excitation  and  ionization.  Energy 
transfer  into  elastic  heating  is  generally  very  small  for  the  range  of  E/N  of 
interest  here,  and  is  therefore  neglected. 

For  reasons  of  both  simplicity  and  generality  it  is  desirable  to  express 
the  various  impact  excitation  rates  in  terms  of  a  parameter  that  renders  them 
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Fig.  9  -  Fractional  Discharge  Energy  Transfer  into  Various  Excitation  Levels 
in  CO,/CO/N?/H,  1/0.25/5.5/0.125 
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Fig,  11  -  Electron  Impact  Vibrational  Excitation  Rate  of  002(100+020) 
for  Various  Gas  Mixtures 
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Fig.  16  -  Electron  Impact  Electronic  Excitation  Rate  for  CO?  and  N7 
in  Various  mixtures  2  2 
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Fig.  19  -  Electron  Impact  Ionization  Rate  of  CO  for  Various 
Gas  Mixtures 
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independent  of  the  gas  mixture.  As  previously  mentioned  by  Elliott  et  al. 

(Ref.  12),  the  average  electron  energy  is  such  a  parameter,  at  least  nearly  so. 
Figures  10  through  12  show  that  vibrational  excitation  rates  for  CC^OIO), 
00^(100  + 020),  and  CC^fOOl),  when  plotted  as  a  function  of  average  electron 
energy,  are  practically  insensitive  to  mixture  variations  over  a  wide  range 
of  different  mixture  compositions.  Figure  13  indicates  that  the  00^(030) 
excitation  is  slightly  faster  for  mixtures  with  He  as  compared  to  those  con¬ 
taining  H£.  Still,  the  difference  is  negligible  in  view  of  other  uncertainties 
of  the  overall  plasma  chemistry  model.  The  ^(1)  excitation  rate  shown  in 
Fig.  14  shows  some  scatter  only  at  the  very  low  electron  energies.  No  clear 
distinction  between  He  and  H^  mixture  is  discernable  here.  Figure  15,  show¬ 
ing  the  CO(l)  vibrational  excitation  rate,  indicates  that  the  presence  of  does 
have  a  noticeable  effect  on  this  process,  the  rate  being  faster  for  mixtures  with 
N£  as  compared  with  a  mixture  without  N^.  All  in  all,  with  regard  to  the  dom¬ 
inant  electron  impact  vibrational  excitation  rates  of  CC>2#  and  CO,  the  use 
of  mixture  independent  excitation  rates  as  a  function  of  electron  energy  should 
be  an  excellent  approximation. 

Unfortunately  the  same  cannot  be  said  about  the  electronic  excitation 
rates  for  CO and  CO  which  arc  shown  in  Figs.l6  and  17.  First  of  all, 
it  should  be  noticed  that  these  rates  are  roughly  at  least  an  order  of  magni¬ 
tude  lower  than  the  vibrational  excitation  rates.  For  average  electron  energies 
below  1.0  eV  they  may  be  several  orders  of  magnitude  lower.  The  rate  of 
N2  excitation  appears  to  be  dependent  on  the  presence  of  He  or  H2  while  the 
relative  amounts  of  other  gases  do  not  seem  to  have  any  noticeable  effect. 

For  electronic  excitation  of  CC^,  the  rate  can  vary  up  to  almost  an  order  of 
magnitude  depending  on  the  gas  mixture,  as  seen  from  Fig.  16.  The  rate  for 
electronic  excitation  of  CO  (Fig.  17)  also  depends  on  the  gas  mixture  com¬ 
position  to  some  degree.  Mixtures  containing  N2  result  in  a  distinctly  lower 
rate  for  this  process  as  compared  to  the  mixture  without  Ng. 

From  the  ionization  cross-section  data  as  supplied  with  the  Boltzmann 
code  (Ref.  8)  it  can  be  inferred  that  ionization  of  He  and  N2  does  not  play  any 
role  except  at  very  high  electron  energies  (u  >  18  eV).  Ionization  of  CO?, 


CO,  and  does  occur  at  rates  which  are  several  orders  of  magnitude  lower 
than  those  for  vibrational  and  electronic  excitation.  Ionization  rates,  unfor¬ 
tunately,  are  very  dependent  on  the  particular  gas  mixture  as  shown  in  Figs. 

18  and  19.  The  oxygen  ionization  rate  in  a  typical  He  mixture  is  shown  in 
Fig. 20. 

The  discussion  of  the  Boltzmann  code  results  permits  the  following 
conclusions.  When  a  new  gas  mixture  is  to  be  analyzed,  Boltzmann  code 
calculations  have  to  be  made  to  determine  the  drift  field  velocity  as  a  func¬ 
tion  of  E/N.  It  is  very  likely  that  no  changes  are  needed  in  the  electron 
impact  vibrational  excitation  rates.  Rate  constants  for  electronic  excitation 
and  ionization  will  most  probably  have  to  be  determined  for  each  individual 
mixture. 

There  is  one  additional  problem  that  must  be  mentioned.  In  the 
Boltzmann  code,  CC>2  (100  +  020)  is  treated  as  a  compound  state  resulting  in 
a  compound  rate  constant  for  impact  vibrational  excitation  of  these  levels. 

Since  the  plasma  kinetics  calculations  consider  CO^(IOO)  and  CC>2(020)  as 
separate  species,  some  prescription  has  to  be  devised  for  splitting  the  com¬ 
pound  rate  constant.  Because  of  Fermi  resonance  these  two  levels  are  closely 
coupled;  in  fact,  many  C07  laser  codes  explicity  assume  equilibrium  for  these 
levels.  Equilibrium  calculations  indicate  that,  in  the  range  of  200  K  <  T  <  400 K, 
the  population  ratio  is  roughly  C02( 100)/C02(020)  =  1/4.  In  the  absence  of 
more  precise  information  we  therefore  assume  that 

k(CO2(100))  =  0.2  k(C02(100  +  020)) 

k(C02(020))  =  0.8k(C02(100  +  020)) 

for  use  in  the  plasma  chemistry  and  laser  calculations.  This  argument  is 
supported  by  comparisons  of  our  analysis  with  experimental  data  to  be  dis¬ 
cussed  in  Section  3. 

*Nmnp/N  =[e-mVT  (n+1)  e  '  n62/T  c  '  P63/T]/Qvlb 
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2.2.3  General  Reaction  Mechanisms 


Aside  from  primary  ionization  reactions  and  secondary  excitation  mech¬ 
anisms  described  in  detail  in  the  preceding  subsections,  EDL  plasma  chemistry 
calculations  require  a  large  number  of  reactions  dealing  with  electron  attach¬ 
ment,  detachment  change  transfer,  and  recombination,  electronic  de-excitation 
and  vibrational  relaxation,  and  finally,  a  sizable  number  of  purely  chemical 
reactions.  The  general  reaction  set  used  for  this  effort  is  an  updated  version 
of  a  set  which  was  previously  documented  in  detail  (Refs.  3  and  13).  A  listing 
of  the  set  currently  used  is  provided  in  the  Appendix. 

2.3  CLOSED  CYCLE  SIMULATION 

One  of  the  requirements  of  the  present  investigation  was  to  further  in¬ 
vestigate  the  response  of  EDL  plasmas  when  subjected  to  recirculation  and 
repetitive  pulsing  in  a  closed  cycle  system.  Several  recent  investigations, 
experimental  as  well  as  theoretical,  have  shown  that  closed  cycle  electric 
discharge  lasers  are  subject  to  the  formation  of  contaminants  which,  in  all 
cases,  tend  to  degrade  laser  performance.  For  example,  electron-beam  in¬ 
teraction  with  materials  may  cause  signigicant  amounts  of  attaching  contami¬ 
nants  to  be  released.  We  assumed  ideal  non- reacting  walls  to  explore  the 
effects  of  discharge  produced  plasma  by-products.  Discharge -induced  con¬ 
taminant  formation  is  a  time  dependent  process  which  greatly  depends  on 
initial  gas  composition,  pressure  level,  discharge  energy  and  pulse  width, 
making  it  extremely  difficult  to  draw  any  conclusions  of  general  validity. 

Kinetics  calculations  are  dependent  on  density,  temperature  and  gas 
composition.  To  predict  the  formation  of  contaminants  it  is  therefore  neces¬ 
sary  to  simulate  the  varying  conditions  that  the  gas  experiences  as  it  traverses 
the  closed  cycle  system.  This  includes  the  consideration  of  catalytic  devices 
designed  to  remove  the  most  important  electrophilic  reaction  products,  there¬ 
by  eliminating  the  problem  of  performance  degradation.  Of  equal  importance 
is  the  consideration  of  the  effect  of  heat  exchangers  designed  to  return  the 
recirculating  gas  mixture  to  the  initial  temperature  after  each  cycle.  Con¬ 
densation  effects,  however,  are  not  considered. 
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It  is  here  assumed  that  the  complete  cycle  can  be  simulated  in  terms 
of  three  major  events,  as  shown  in  Fig.  21.  The  first  leg  of  the  triangular 
cycle  in  the  p-V  diagram  represents  the  discharge  itself.  The  discharge  is 
followed  by  an  adiabatic  expansion.  It  is  assumed  that  after  the  expansion 
the  gas  is  again  at  the  initial  pressure.  The  temperature  at  that  point  is 
generally  considerably  higher  than  the  pre-pulse  temperature,  and  even 
higher  than  the  post-pulse  temperature  due  to  collisional  relaxation  of  the 
energy  stored  in  the  vibrational  levels  during  the  discharge  pulse.  While  the 
time  frame  for  the  discharge  is  given  by  the  discharge  pulse  duration,  the 
duration  of  the  adiabatic  expansion  Ls  estimated  by  calculating  the  time  inter¬ 
val  which  it  would  take  for  a  pressure  pulse  to  clear  the  particular  cavity 
being  considered.  The  final  leg  of  the  triangular  cycle  simulates  the  return 
of  the  mixture  to  the  initial  conditions  (pressure  and  temperature).  This  leg 
represents  a  relaxation  calculation  at  constant  pressure  and  specified  distri¬ 
bution  of  temperature  versus  time.  Typically,  the  time  frame  for  this  leg 
represents  the  major  portion  of  the  gas  transit  time  through  one  complete 
cycle.  The  result  of  the  relaxation  calculations  through  this  triangular  path 
is  that  the  gas  is  now  at  the  initial  (pre-pulse)  pressure  and  temperature, 
but  has  a  composition  which  is  generally  different  from  that  at  the  beginning 
of  the  cycle. 

Removal  of  contaminants  via  catalyst  action  is  currently  considered  at 
this  point.  The  effect  of  a  drying  bed  included  to  remove  water  vapor  being 
generated  in  certain  gas  mixtures  can  be  modeled  in  the  same  fashion.  Know¬ 
ing  the  gas  composition  at  the  end  of  the  cycle,  the  effect  of  the  catalytic  reactor 
and/or  drying  bed  is  computed  by  adjusting  all  concentrations  according  to 

F!  =  (1  -  rji )  F.  (32) 

where  represents  the  catalyst/drying  bed  removal  efficiency  for  the  ith 
species,  and  F^  its  mole/ mass  ratio  upstream  of  the  catalyst/drying  bed.  The 
cycle  is  then  repeated  starting  with  the  adjusted  composition  F|  at  the  initial 
temperature  and  pressure. 
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3.  CODE  VALIDATION 


3.1  INTRODUCTION 

It  is  obviously  very  important  to  gain  some  confidence  in  the  validity  of 
the  analysis  described  in  Section  2.  Comparison  of  the  analytical  model  with 
experimental  data  will  give  a  direct  assessment  of  the  model  as  a  predictive 
tool  for  designing  closed  cycle  EDL  systems  for  optimum  performance.  Al¬ 
though  earlier  versions  of  the  analytical  model  (Ref.  14)  have  shown  good  over¬ 
all  agreement  with  experimental  data  (small  signal  gain,  discharge  power 
loading,  laser  efficiency),  numerous  changes  and  refinements  that  have  been 
incorporated  into  the  model  in  the  course  of  the  present  effort  make  it  neces¬ 
sary  that  further  comparisons  between  theory  and  experiment  be  made.  While 
the  present  effort  was  originally  undertaken  with  the  goal  of  defining  optimum 
operating  conditions  for  two  closed  cycle  EDL  systems  presently  under  devel¬ 
opment,  delays  in  the  completion  of  these  devices  have  made  it  impossible  to 
obtain  any  experimental  data  from  these  devices. 

Fortunately  some  other  data  have  been  made  available  to  us.  The  first 
case  to  be  discussed  in  this  section  concerns  a  small  scale  closed  system 
operated  at  the  MIT  Lincoln  Laboratories  with  the  purpose  of  investigating 
the  effects  of  discharge  generated  contamination  on  performance  as  a  function 
of  time  (Ref.  15). 

The  second  case  concerns  a  large  open  cycle  device  (ABEL).  Even 
though  discharge  generated  contamination  Is  of  lesser  importance  in  an  open 
cycle  device,  it  afforded  us  the  opportunity  to  compare  model  results  with 
projected  performance  data  (Ref.  16). 
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3.2  SMALL  SCALE  E-BEAM  DEVICE 

Reportedly,  the  closed  system  small  scale  E-beam  device  tested  at  the 
Lincoln  Laboratories  was  designed  with  the  aim  to  minimize  contamination 
of  the  laser  gas  from  the  device  hardware  (outgassing,  etc.).  The  nominal 
operating  conditions  are  summarized  in  Table  2,  and  Fig.  22  shows  a  cross 
section  of  the  hardware  used.  Experimental  results  were  provided  in  the 
form  of  oscilloscope  traces  for  the  sustainer  current  and  the  small  signal 
gain  of  the  P(20)  transition.  As  no  heat  exchanger  was  used,  heat  conduction 


Table  2 

CLOSED  SYSTEM,  SMALL  SCALE  E-BEAM  DEVICE  DATA 
(Lincoln  Laboratory,  Ref.  15) 


Gas  Composition 
P  (Initial) 

T  (Initial) 

E-Beam  Current 
E-Beam  Energy 
Sustainer  Voltage 
Electrode  Gap 
Cross  Section 
Gain  Length 
Pulse  Duration 
Pulse  Frequency 
Post-Pulse  Expansion 
Cycle  Time 


C02/C0/N2/He/02 
605  torr 
298  K 

1.6  mA  (Effective) 

1 65  keV  (Nominal 
8.1  kV 
3  cm 

See  Fig.  22 
15  cm 
300  /isec 
1  Hz 

300  ji sec 
1.0  sec 


1/1/4/6/0.1 


to  the  water-cooled  cavity  walls  was  relied  upon  to  reduce  the  post-pulse 
temperature  to  room  temperature  between  pulses.  Simplified  heat  conduction 
calculations  indicated  that  this  method  should  indeed  suffice. 

In  any  attempt  to  model  an  experiment  which  is  not  tailored  to  the  analy¬ 
sis  certain  assumptions  will  have  to  be  made  before  actual  calculations  can  be 
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performed.  For  example,  examination  of  the  cavity  cross  section  (Fig.  22) 
shows  that  discharge  conditions  in  the  cavity  are  not  likely  to  be  spatially 
uniform  as  assumed  by  the  theoretical  model.  In  this  case  we  have  chosen 
to  model  the  center  section  of  the  cavity,  i.e.,  a  region  which  has  the  same 
width  as  the  foil  window,  and  which  includes  the  location  of  the  probe  laser 
beam.  This  is  equivalent  to  assuming  a  collimated  electron  beam,  or  neg¬ 
lecting  electron  beam  spreading  in  the  cavity  due  to  whatever  reason. 

A  quantity  that  requires  further  examination  is  the  electron  beam  energy. 
Beam  electrons  experience  losses  due  to  interaction  with  the  foil  and  due  to 
collisions  in  the  cavity.  Some  of  these  losses  are  compensated  for  by  accelera¬ 
tion  in  the  electric  field  of  the  discharge.  Using  the  data  by  Berger  and  Seltzer 
(Ref.  7)  the  energy  loss  in  the  foil  under  the  given  conditions  is  estimated  to 
be  about  17  keV.  The  average  collisional  loss  in  the  cavity  amounts  to  another 
2  keV,  while  the  discharge  will  provide  an  average  gain  of  about  4  keV.  This 
results  in  a  net  loss  of  about  15  keV,  and  therefore  yields  an  effective  electron 
beam  energy  of  165  -  15  =  150  keV. 

Finally,  in  order  to  perform  the  plasma  chemistry  calculations,  the 
electron  impact  excitation  rates  have  to  be  established  using  the  Boltzmann 
code.  These  calculations  and  their  results  have  already  been  discussed  in 
Section  2. 

All  in  all,  only  relatively  minor  corrections  were  needed  in  the  model 
to  arrive  at  the  comparison  between  experimental  and  theoretical  results 
shown  in  Figs.  23  and  24.  The  major  problem  encountered  was  the  uncertainty 
in  the  dominant  rate  of  electron  attachment  to  O,.  While  this  rate  was  un- 
important  in  previously  computed  cases,  it  turned  out  to  be  of  great  importance 
in  this  particular  gas  mixture  (see  Table  2).  Reducing  the  magnitude  of  the 
rate  constant  by  a  factor  of  two  (which  was  well  within  the  range  of  uncertainty) 
and  giving  it  an  inverse  temperature  dependence  (as  suggested  by  theoretical 
considerations)  the  temporal  distributions  of  small  signal  gain  and  sustainer 
current  as  calculated  are  in  satisfactory  agreement  with  the  experimental  data 
as  shown  in  Figs.  23  and  24,  respectively.  The  reason  for  the  discrepancy  in  the 
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Discharge  Current  (A) 
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Fig.  24  -  Calculated  Temporal  Distribution  of  Discharge  Current 
(5th  Pulse) 
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tail  of  the  gain  distribution  is  not  known  at  this  time.  Calculations  with  slightly 
increased  electron  beam  current  have  shown  that  this  could  be  remedied  but 
only  at  the  expense  of  a  higher  peak  gain  and  a  lesser  agreement  between 
theoretical  and  measured  discharge  current  distribution. 

As  mentioned  above,  the  purpose  of  this  experiment  was  to  investigate 
discharge  generated  contamination  and  its  effect  on  discharge  power  loading. 
The  analysis  predicts  the  major  by-products  to  be  various  nitrogen  oxides  as 
shown  in  Fig.  25.  For  the  dry  gas  mixture  the  effect  of  these  nitrogen  oxides 
on  the  discharge  energy  loading  is  minimal  as  can  be  inferred  from  Fig.  26. 

Computer  run  time  for  the  dry  gas  case  with  ten  cycles  was  roughly 
6  minutes  on  the  CDC  Cyber  176.  This  shows  clearly  that  it  would  be  pro¬ 
hibitive  to  continue  these  calculations  for  something  like  2000  cycles.  It  is 
therefore  of  interest  to  develop  an  estimate  of  the  performance  decay  rate  as 
a  function  of  time,  and  to  compare  the  results  with  experimental  observations. 
As  will  be  seen,  it  is  then  also  of  importance  to  know  that  the  "gas  utilization 
factor"  in  this  experiment  was  about  1/ 1 5  (plasma  volume  to  total  system  gas 
volume).  Assuming  an  exponential  decay  for  the  discharge  power  loading  of 
the  form 


E. 


in 


A  • 


e 


-an 


(33) 


where  n  denotes  the  number  of  cycles  (or  pulses  in  this  case),  the  time  con¬ 
stant  or  decay  rate,  a  can  be  determined  from  the  results  shown  in  Fig.  26. 

_  3 

This  yields  a  value  of  a  =  1.31  •  10  Considering  the  gas  utilization  factor 
of  g  =  1/ 15 ,  the  actual  decay  rate  would  be  approximately  a  •  g  =  8.7  •  10 
With  this  value,  we  find  from  Eq.  (33)  for  n  =  2000 


e-a  •  g  •  2000 


0.84 


(34) 
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Pulse  Number 


Normalized  Discharge  Power  Loading  as  a  Function 
of  Pulse  Number  (C02/C0/N2/He/C>2  =  1/1/4/6/0.1) 


which  indicates  an  overall  decay  of  16%  after  2000  pulses.  This  is  in  excel¬ 
lent  agreement  with  the  experiment  which  showed  a  decay  of  10  to  20%  in  (he 
gain  after  2000  pulses  (Ref.  15). 

Performance  deteriorates  drastically  when  water  vapor  is  added  to  the 
mixture.  The  effect  of  30  ppm  H^O  on  the  discharge  energy  loading  is  also 
shown  in  Fig.  26.  The  simultaneous  presence  of  H^O  (or  H^)  and  certain 
nitrogen  oxides  (see  Fig.  27)  generally  leads  to  the  formation  of  HNO^,  NHO^ 
and  HO^NO^  (see  Fig.  28)  all  of  which  have  large  cross  sections  for  electron 
attachment.  Although  the  HNOx  concentrations  shown  in  Fig.  28  are  of  the 
same  magnitude  as  the  NOx  concentrations  shown  in  Fig.  27,  electron  attach¬ 
ment  cross  sections  for  the  HNO^  are  several  orders  of  magnitude  larger  than 
those  for  NOx>  thereby  accounting  for  the  decay  in  discharge  current  and  small 
signal  gain  (Fig.  29)  when  moisture  is  present  in  the  discharge  gas  mixture. 

3.3  LARGE  OPEN  CYCLE  DEVICE  (ABEL) 

ABEL  is  a  cold  flow,  open  cycle,  pulsed  electron  beam  sustained  electric 
discharge  laser.  Its  principal  purpose  is  to  develop  and  prove  high  pulse  rate, 
high  power  electric  discharge  laser  technology.  A  schematic  of  the  cavity 
configuration  is  shown  in  Fig.  30.  Cavity  characteristics  and  typical  operating 
conditions  as  well  as  projected  performance  data  are  listed  in  Table  3.  As 
usual,  the  flow  direction,  the  direction  of  the  electric  field  and  the  optical  axis 
are  mutually  perpendicular. 

Electron  impact  excitation  rates  for  this  case  were  discussed  in  Section  2. 
The  only  additional  assumption  used  in  the  calculations  for  this  device  concerns 
the  electric  field  distribution  which  is  assumed  to  have  a  linear  rise  from  one 
tenth  of  nominal  value  to  full  nominal  value  over  the  first  two  microseconds  of 
the  discharge  pulse. 

Some  illustrative  results  are  presented  in  Figs.  31  through  33,  showing 
the  small  signal  gain  for  the  P(18)  transition,  temperature  distributions  with 
and  without  lasing,  and  the  secondary  electron  density,  respectively. 
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Fig,  3  1  -  P(  18)  Small  Signal  Gain  Temporal  Distribution 


CO,/N,/H,  =  1/3/0. 1  (ABEL) 
E/N  =  3.7  •  10"  1  Vcrn 

•>EB  =  60  A 


With  Lasing 


J _ I _ I _ I _ I _ L 

4  8  12  16  20  24 

Time,  t  (/Xsec) 


Fig.  33  -  Secondary  Electron  Density  Temporal  Distribution 


The  most  important  predictions  are,  however,  those  of  output  flux  in¬ 
tensity  and  discharge  power  loading  and  their  temporal  distribution  during 
pulsed  laser  operation.  The  P(18)  output  flux  distribution  shown  in  Fig.  34 
represents  an  output  power  of  0.15  kj/lit.  The  corresponding  discharge 
power  loading  distribution  shown  in  Fig.  35  represents  a  discharge  energy 
input  of  1.09  kj/lit  with  lasing,  and  0.98  kj/lit  without  lasing.  The  result¬ 
ing  efficiency  therefore  is  about  14%.  While  the  output  power  is  in 
excellent  agreement  with  the  projected  performance,  the  discharge  power 
loading  appears  to  be  about  20  to  25%  high,  resulting  in  an  efficiency  which 
is  probably  on  the  low  side.  The  dependence  of  all  three  quantities  on  the 
effective  electron  beam  current  is  shown  in  Fig.  36.  It  can  be  seen  that  at 
a  current  of  40  Amp  the  discharge  power  loading  becomes  about  0.85  kj/lit 
and  the  output  pulse  energy  decreases  slightly  to  about  0.14  kj/lit  with  a 
resulting  efficiency  of  16%.  Therefore,  making  allowance  for  the  uncer¬ 
tainty  in  the  given  effective  electron  beam  current  it  appears  that  the  present 
analysis  is  in  excellent  agreement  with  the  projected  performance  of  the  de¬ 
vice  under  study. 

Finally,  it  should  be  remembered  that  these  calculations  were  per¬ 
formed  for  given  pulse  durations  (see  Table  3).  It  is  clear  from  Fig.  34  that 
the  pulse  stretches  as  the  electron  beam  current  (and  the  secondary  electron 
density  as  a  consequence)  is  decreased.  Figure  34  shows  that  for  E-beam 
currents  of  less  than  100A  the  output  pulse  is  terminated  by  the  E-beam  cut¬ 
off.  Conclusions  drawn  from  Fig.  36  must  keep  this  in  mind.  In  future  calcu¬ 
lations  the  pulse  duration  should  be  considered  a  variable  to  be  determined 
as  a  result  of  the  analysis  rather  than  specified  a  priori. 
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CO,/N,/H,  =  1/3/0. 1  (ABEL) 
E/N  =  3.7  *10" 16  Vcm2 


4.  ANALYSIS  OF  TEST  CASES 

4.1  AFWL  PULSAR  MODELING 

The  PULSAR  device  is  a  closed  cycle,  electron-beam  sustained  electric 
discharge  laser  presently  under  development  at  AFWL  for  the  purpose  of 
studying  discharge  stability,  medium  homogeneity,  acoustics,  energy  loading 
and  general  discharge  physics  under  closed  cycle  operation.  Also  planned  is 
the  investigation  of  fast  startup  transients.  The  purpose  of  this  section  is  to 
describe  a  modeling  effort  specifically  aimed  at  investigating  the  effects  of 
discharge  generated  contamination  of  PULSAR  closed  cycle  operation. 

4.1.1  Open  Cycle  Performance 

Before  directing  our  attention  to  closed  cycle  operation  it  is  of  interest 
to  establish  a  reference  frame  in  terms  of  open  cycle  operation.  Table  4 
lists  some  basic  given  characteristics  for  this  device,  such  as  nominal  gas 
mixture,  initial  pressure  and  temperature,  maximum  discharge  energy  loading 
and  basic  cavity  geometric  data.  The  pulse  duration  times  were  assumed 
based  on  given  data  for  the  Army  Small  Circulator  to  be  discussed  later.  Initial 
calculations  also  assumed  a  discharge  potential  of  25  kV  and  an  effective  elec¬ 
tron  beam  current  density  of  100  mA/cm^,  values  which  were  specified  for  the 
Army  Small  Circulator.  The  indicated  data  yield  an  electric  field  to  density 
ratio  of  E/N  =  1.36  •  10  ^  Vcm^,  a  value  which  according  to  Fig.  8  is  much 
too  low  to  achieve  efficient  fractional  energy  transfer  into  the  relevant  vibra¬ 
tional  levels.  In  fact,  Fig.  8  indicates  an  optimum  value  of  E/N  %  3.5  •  10  ^ 

— »4.0  ♦  10  ^  Vcm^.  To  stay  below  a  value  at  which  measurable  discharge 
energy  is  consumed  in  the  processes  of  electronic  excitation  and  ionization  we 
chose  a  discharge  potential  of  65  kV  which  yields  a  value  of  E/N  =  3.54  •  10-1^ 
Vcm  ,  very  close  to  the  optimum  as  indicated  by  the  Boltzmann  code  results. 

In  reality,  of  course,  the  requirement  to  avoid  arcing  will  impose  the  upper 
limit  to  this  value.  Results  of  discharge  laser  performance  calculations 
(arbitrarily  assuming  a  50%  Fabry-Perot  resonator  output  coupling)  as  a 
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Table  4 

AFWL  PULSAR  DEVICE  DATA 
(Ref.  17) 


Gas  Composition  (Nominal) 

COz/N2/H2  =  1/3/0.08 

p  (Initial) 

760  torr 

T  (Initial) 

200  K 

Discharge  Energy  Loading 

800  j/lit  (max) 

Electrode  Gap 

5  cm 

Cavity  Volume 

5  x  5  x  25  cm 

Gain  Length 

25  cm 

Pulse  Duration 

E-Beam 

8  fi  sec 

Discharge 

10  fisec 

function  of  E/N  are  shown  in  Fig.  37.  It  becomes  clear  then,  that  operation 
at  optimum  E/N  with  an  effective  electron  beam  current  of  100  mA/crn^  for 
the  assumed  pulse  duration  given  in  Table  4  far  exceeds  the  energy  loading 
limit  of  the  PULSAR  hardware.  Since  discharge  energy  loading,  for  recom¬ 
bination  dominated  plasmas,  is  roughly  proportional  to  the  square  root  of  the 
effective  electron  beam  current,  decreasing  the  effective  electron  beam  cur¬ 
rent  while  maintaining  the  optimum  electric  field  should  lower  the  discharge 
energy  loading  to  an  acceptable  value.  The  results  of  these  calculations  are 
shown  in  Fig.  38,  indicating  an  upper  limit  of  around  30  mA/cm^  for  the 
effective  electron  beam  current  density.  The  output  energy  and  the  efficiency 
vary  quite  strongly  as  a  function  of  the  effective  electron  beam  current,  as 
indicated  in  Fig.  39.  It  should  be  noted  that  lowering  the  effective  E-beam 
current  density  from  100  mA/cm  to  30  mA/cm  also  lowered  the  output 
energy  by  approximately  25%,  but  increased  the  efficiency  by  roughly 
75%.  To  illustrate  further  the  effect  of  E-beam  current  density  on  the 
discharge,  secondary  electron  densities  and  gas  temperatures  are  shown  as 
a  function  of  time  during  the  pulse  in  Figs.  40  and  41,  respectively.  Sum¬ 
marizing  the  discussion  up  to  this  point,  we  have,  for  the  conditions  given  in 
Table  4,  determined  the  optimum  electric  field  strength  from  the  results  of 
the  Boltzmann  code.  For  this  value,  variations  in  the  effective  electron 
beam  current  density  have  shown  that  a  current  density  of  30  mA/cm^  is 
consistent  with  the  discharge  energy  loading  capability  of  the  PULSAR  dis¬ 
charge  hardware.  We  are  now  in  a  position  to  investigate  the  response  of 
the  system  to  variations  in  gas  mixture  composition. 

It  is  of  particular  interest  to  study  variations  in  the  concentration 
of  the  nominal  CO^/N^/H^  =  1/3/0.08  mixture.  Since  is  needed  to  provide 
V-T  relaxation  of  the  CO^  bending  mode,  thus  providing  a  drain  for  the  lower 
laser  level  population,  the  presence  of  also  is  responsible  in  this  system 
for  the  formation  of  electrophilic  by-products  such  as  HNO^,  HNO^  and 
HO^NO^.  The  question  arises;  can  a  reduction  in  concentration  minimize 
the  contaminant  formation  without  adversely  effecting  putput  power  and  lasing 
efficiency? 
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Fig.  41  -  Gas  Temperature  as  Function  of  Time  for  Varying 
E-Beam  Current  Density  (Open  Cycle,  with  Lasing) 
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Again,  it  is  illustrative  to  study  this  effect  for  open  cycle  (or  single  pulse) 
operation.  It  is  found  that  the  gas  mixture  indeed  shows  optimum  performance 
for  0.08  parts  as  shown  in  Fig.  42.  Both  output  energy  and  efficiency  show 
a  broad,  flat  optimum  with  relatively  little  change  between  0.05  and  0.16  parts 
H^.  Lowering  the  content  below  0.05  will  cause  a  strong  decrease  in  both 
output  energy  and  efficiency.  From  Fig.  43,  which  shows  the  corresponding  dis¬ 
charge  current  density  distributions  as  a  function  of  time  and  content,  it  is 
concluded  that  discharge  energy  loading  is  not  much  affected  by  changes  in  the 
concentration  for  content  varying  between  0.02  and  0.16.  For  0.005  the 
energy  input  is  about  6.8%  higher  than  for  0.08.  Output  intensity  distributions  of 
the  P(l6)  transition  as  a  function  of  content  are  shown  in  Fig.  44.  Figure  45 
shows  some  typical  discharge  current  densities  without  the  effect  of  lasing.  It  is 
noted  that  a  slight  increase  in  current  density  is  obtained  as  the  content  is 
lowered  from  0.08  to  0.02.  It  should  be  noted  that  it  is  difficult  to  change  any 
single  parameter  without  affecting  others  and  this  may  provide  the  explanation 
here.  As  the  content  is  reduced  from  0.08  to  0.02  at  constant  pressure,  the 

density  increases  by  1.4%  while  the  discharge  power  input  increases  by  2.64%. 

2 

Since,  to  a  first  approximation,  the  electron  density  is  proportional  to  p  ,  squaring 
the  density  yields  A p  ~  2.8%, 'hence  the  correspondingly  higher  current  density. 
Concluding  the  open  cycle  analysis.  Fig.  46  shows  small  signal  gain  distributions 
for  varying  concentrations.  While  the  amount  of  does  not  greatly  affect 
the  small  signal  gain,  it  does  noticeably  affect  the  open  cycle  lasing  perform¬ 
ance  as  already  shown  in  Figs.  42  and  44. 

While  we  conclude  from  Fig.  43  that  discharge  energy  loading  is  nearly 
independent  of  concentration  (for  the  range  from  0.02  to  0.16  parts  H^),  it 
is  noted  that  over  the  same  range  lasing  causes  the  discharge  energy  loading 
to  be  between  4  and  6%  higher  than  without  lasing  as  shown  in  Fig.  47.  For 
open  cycle  operation  the  formation  of  contaminants  in  the  discharge  should  be 
important  only  if  it  affects  the  electron  density  during  the  actual  pulse.  Figure 
48  shows  that  lasing  has  a  negligible  effect  on  open  cycle  plasma  by-product 
formation  except  for  NO^.  The  likely  principal  cause  for  the  increase  is  dis¬ 
charge  energy  loading  with  lasing  is  the  observed  faster  deactivation  of  high 
vibrationally  excited  levels  (particularly  N^l),  see  Fig. 49).  This  leads  to 
a  higher  gas  temperature  which  effects  the  rates  such  that  the  electron  con¬ 
centration  is  increased  as  compared  to  the  non-lasing  case. 


P(  1 6 )  Output  Energy  (J/lit  atm) 


Fig.  42  -  P(16)  Output  Energy  and  Efficiency  as  Function 
of  H£  Content  (P(16),  Open  Cycle) 


Efficiency  (%) 


Fig.  43  -  Discharge  Current  as  Function  of  Time  and 
Concentration  (Open  Cycle,  with  Lasing) 
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Fig .  46  -  Small  Signal  Gain,  P(16)  as  F unction  of  Time  and  H£ 
Concentration  (Open  Cycle) 
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Fig.  47  -  Effect  of  Lasing  on  Discharge  Energy  Loading 
(Open  Cycle) 
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Fig.  48  -  Effect  of  Lasing  on  Contaminant  Generation  During 
Discharge  Pulse  (Open  Cycle)  for  CO?/N,/H,  = 
1/3/0.08  (PULSAR)  tie. 
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Effect  of  Lasing  on  Dominant  Species  (Open  Cycle) 


4.1.2  Closed  Cycle  Performance 

The  simulation  of  closed  cycle  operation  is  performed  by  continuing  the 
plasma  chemistry  calculations  beyond  the  discharge  pulse  around  the  recircu¬ 
lation  cycle  as  described  in  Section  2.3.  Based  on  typical  post- pulse  gas  tem¬ 
peratures  and  the  given  cavity  geometry  the  time  for  the  adiabatic  expansion 
was  estimated  to  be  about  100  /isec,  and  this  value  was  used  throughout  the 
closed  cycle  performance  modeling  of  this  case. 

Initial  calculations  used  a  time  frame  of  At  =  5  msec  with  a  linear  de¬ 
crease  of  the  temperature  to  simulate  the  effect  of  a  heat  exchanger  in  the 
remainder  of  the  cycle.  Since  the  results  showed  still  very  strong  species 
gradients  at  the  end  of  this  assumed  cycle,  the  time  frame  was  arbitrarily 
extended  to  At  =  1.0  sec.  The  effect  of  different  cycle  times  on  normalized 
discharge  energy  loading  for  five  cycles  is  shown  in  Fig.  50.  Not  only  had  all 
species  gradients  become  vanishingly  small  for  At  =  1.0  sec,  the  power  loading 
decay  was  substantially  less  than  for  the  shorter  time  frame.  These  results 
shew  that  it  is  rather  important  to  employ  a  realistic  history  of  temperature 
versus  time  in  the  simulation  of  the  recirculation  loop. 

Realistic  gasdynamic  conditions  around  the  PULSAR  recirculation  loop 
were  obtained  from  AFWL  (Ref.  17).  Table  5  lists  their  estimates  for 
velocity  and  temperature  as  a  function  of  distance  around  the  tunnel.  Static 
pressure  values  were  found  to  show  only  small  variations  and  were  therefore 
ignored.  The  main  interest  here  is  the  temperature,  and  its  effect  on  density. 
Logarithmic  differentiation  of  the  equation  of  state  shows  that  temperature 
variations  have  a  far  more  pronounced  effect  on  density  than  small  pressure 
variations.  To  obtain  a  better  picture  of  the  recirculation  loop,  velocity  and 
temperature  from  Table  5  are  plotted  versus  distance  in  Figs.  51  and  52, 
respectively.  To  establish  the  temperature  versus  time  history,  the  reciprocal 
velocity  is  integrated  versus  flow  distance  assuming  a  linear  variation  between 
given  data  points.  The  results  are  tabulated  in  Table  6  and  displayed  in  Fig.  53, 
along  with  a  typical  temperature  distribution  during  the  discharge,  adiabatic 
expansion  and  early  part  of  the  cycle  up  to  the  first  heat  exchanger.  From  this 
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COz/N2/H2  =  1/3/0.08  (PULSAR) 
E/N  =  3.54  *10' 16  Vcm2 
jEB  =  30  mA/cm2 
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Tabic  5 

AFWL  PULSAR  DEVICE 
RECIRCULATION  LOOP  DATA 


Station 

X 

(cm) 

u 

(m/s) 

B 

Remarks 

0.0 

24.4 

208 

Cavity  Entrance 

i 

5.0 

64.5 

541 

Cavity  Exit 

2 

35.5 

64.5 

541 

End,  Hot  Duct 

3 

85.5 

25.0* 

541 

End,  Diffuser 

4 

100.0 

25.0 

541 

Heat  Exchanger  1 

5 

130.5 

25.0 

284 

End,  Heat  Exch.  1 

6 

145.0 

24.0 

284 

End,  Duct 

7 

178.0 

65.0 

284 

End,  Nozzle 

8 

mms 

65.0 

284 

End,  Duct 

9 

H 

74.5 

350 

End,  Blower 

10 

■9 

74.5 

350 

End,  Duct 

11 

337.0 

15.0 

350 

End,  Diffuser 

12 

437.0 

18.0 

350 

End,  Duct 

13 

472.0 

15.0 

208 

End,  Heat  Exch.  2 

14 

487.0 

15.0 

208 

End,  Duct 

15 

521.0 

24.4 

208 

End,  Nozzle 

16 

560.0 

24.4 

208 

End,  Duct 

Value  estimated  based  on  diffuser  area  ratio. 
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Table  6 
RESULTS  OF  At  =  / 


Station 


* 

(sec) 

0.0 

1.21 

io-3 

5.94 

10*3 

1.79 

10'2 

2.39 

io*2 

3.59 

10'2 

4.19 

io'2 

4.98 

10‘2 

5.32 

io-2 

6.18 

io*2 

6.39 

io-2 

8.03 

10‘2 

1.47 

io-1 

1.70 

io'1 

1.80 

io"1 

1.98 

io'1 

2.14 

io'1 

Temperature  to  be  computed  by  code 


point  onward  the  kinetics  calculations  simulate  diabatic  flow  with  the  tempera¬ 
ture  prescribed  as  shown  in  Fig.  53.  It  is  also  recognized  that  the  gas  tem¬ 
perature  in  the  hot  duct  and  at  the  entrance  of  the  first  heat  exchanger  as 
calculated  here  is  in  close  agreement  with  that  resulting  from  the  AFWL  gas- 
dynamic  cycle  analysis.  Table  6  indicates  a  total  cycle  time  of  about  0.22  sec. 
The  discharge  energy  loading  as  a  function  of  recirculation  using  die  computed 
cycle  time  of  At  *  0.22  sec  is  also  shown  in  Fig.  50,  falling  in  between  the  re¬ 
sults  for  the  previously  assumed  cycle  times. 

Having  constructed  a  reasonably  accurate  model  for  closed  cycle  plasma 
chemistry  analysis  we  can  now  attempt  to  determine  die  effect  of  variations  in 
content  of  the  basic  CC^/N^/H^  =  1/3/0.08  mixture  on  closed  cycle  laser 
performance.  Calculations  have  been  carried  out  for  up  to  15  cycles,  with  and 
without  lasing.  A  typical  run  time  on  the  AFWL/CDC  Cyber  176  for  one  case 
is  about  11  minutes  with  lasing,  and  slightly  less  without  lasing.  Results  for 
the  discharge  energy  loading  as  a  function  of  recirculation  and  content  are 
shown  in  Fig.  54.  These  results  show  that  discharge  energy  loading  is  ob¬ 
viously  higher  with  lasing  as  compared  to  die  case  without  lasing,  just  as  for 
open  cycle  operation.  Note  that  the  rate  of  decay  increases  as  the  content 
of  the  mixture  is  decreased. 

The  effect  of  reduced  hydrogen  content  and  recirculation  on  the  laser 
power  output  is  shown  in  Figs.  55  and  56.  As  can  be  seen,  reducing  the 
hydrogen  content  causes  a  decrease  in  power  output,  and  this  response  is  the 
same  as  for  open  cycle  operation.  Figure  56  also  indicates  that,  with  recircu¬ 
lation  it  is  even  more  important  than  for  open  cycle  operation  to  keep  die 
hydrogen  concentration  close  to  the  optimum,  since  the  rate  of  decay  is  in¬ 
creased  as  the  hydrogen  concentration  is  reduced. 

Figure  57  illustrates  the  build  up  of  the  major  by-products  during  re¬ 
circulation  with  lasing  both  for  the  nominal  mixture  and  for  one  with  reduced 
hydrogen  content.  The  principal  result  is  that  reducing  the  hydrogen  content 
will  apparently  lead  to  increased  concentrations  of  NO,  HNO^  and  HNO^ 
thereby  increasing  the  rate  of  electron  attachment,  which  then  results  in 
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Fig.  54  -  Discharge  Energy  l  oading  as  Function  of  Recirculation 
and  Content  With  and  Without  Lasing 
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decreased  output  power.  It  is  interesting  to  note  that  build  up  of  H^O,  N20 
and  HNC>2  shows  very  little  dependence  on  lasing  in  contrast  to  NO  and  HNO^. 
For  the  case  without  lasing,  shown  in  Fig.  58,  the  C02/N2/H2  =  1/3/0.08 
mixture  shows  almost  three  times  the  amount  of  NO  and  about  four  times  the 
amount  of  HNO^  after  15  pulses  as  compared  to  the  case  with  lasing  (Fig.  57). 

The  foregoing  results  have  been  obtained  without  considering  catalytic 
removal  of  contaminants  nor  have  plasma  chemistry  effects  in  the  cathode 
sheath  been  incorpoated.  While  it  is  now  generally  recognized  that  long  term 
closed  cycle  operation  may  require  the  use  of  catalytic  reactors  in  the  cycle, 
the  analysis  presented  here  allows  some  quantitative  observations  to  be  made 
concerning  this  subject.  As  has  been  mentioned  already,  to  extend  the  plasma 
chemistry  calculations  beyond  15  to  20  cycles  becomes  extremely  expensive 
in  terms  of  computer  run  time.  We  shall  therefore  try  to  arrive  at  an  esti¬ 
mate  of  the  long  term  behavior  of  the  device. 


Assuming  again  (see  Section  3.2)  an  exponential  decay  of  performance, 
we  can  determine  a  rate  of  decay  for  the  laser  output  power  as  a  function  of 
cycle  number.  Using  the  curve  for  the  1/3/0.08  mixture  shown  in  Fig.  55, 
this  results  in 


Eout  =  134  79  6 


-1.19  •  10 


-3 


n 


(35) 


which  gives  an  excellent  fit  for  n  =  5  to  15,  n  denoting  the  cycle  number. 
Equation  (35)  yields  a  "half-life  time"  of  582  pulses  with  the  usual  meaning 
that  after  each  582  cycles  the  laser  output  power  is  reduced  by  50%.  This 
calculation  does  not  consider  a  gas  utilization  factor,  and  is  used  only  to 
illustrate  the  fact  that  without  using  a  catalytic  reactor  to  remove  contam¬ 
inants  such  as  NO,  N20,  HN02  and  HNO^  it  will  be  difficult,  if  not  impossible, 
to  maintain  performance  for  any  length  of  time. 

Since  catalytic  reactors  presently  considered  usually  require  relatively 
high  temperatures  to  operate  efficiently,  the  natural  place  for  such  a  device 
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End  of  Cycle  Number 


Contaminant  Buildup  During  Closed  Cycle 
Operation  Without  Lasing  (PULSAR) 


in  the  recirculation  loop  would  be  just  downstream  of  the  discharge  where  the 
gas  reaches  its  highest  temperature.  Figure  59  shows  temperature  distribu¬ 
tions  around  the  cycle,  with  and  without  lasing,  for  1/3/0.08  and  1/3/0.02 
mixtures.  Also  indicated  is  the  "temporal  location"  of  the  hot  duct  down¬ 
stream  of  the  cavity.  Figure  60  shows  some  typical  distributions  of  major 
contaminants,  as  a  function  of  time  around  the  cycle.  Note  that  HN02  and 
HNOj  are  generally  destroyed  in  the  discharge  and  re-formed  during  die 
post-pulse  relaxation  period.  NO,  N20  and  HzO  maintain  an  almost  constant 
distribution,  whereas  OH  and  N02  are  mainly  generated  in  the  discharge,  and 
then  serve  as  principal  reactants,  in  combination  with  available  NO  and  H20, 
for  the  post-pulse  regeneration  of  HN02  and  HNO^.  This  is  the  environment 
in  which  the  catalytic  reactor  will  have  to  function. 

The  analysis  of,  and  predictions  for,  the  PULSAR  device  will  be  con¬ 
cluded  with  an  assessment  of  the  production  of  electrophilic  species  in  the 
cathode  fall  region.  The  C02/N2/H2  =  1/3/0.08  mixture  without  lasing  will 
serve  as  an  example  (see  Fig.  60).  From  Fig.  54  it  is  found  that  the  discharge 
energy  loading  for  this  case  is  typically  about  700  j/lit  (1  pulse)  which,  using 
the  given  cavity  geometry,  translates  into  an  average  discharge  current  den¬ 
sity  ofJr,  =  5.4  A/cm2  during  the  pulse.  Using  a  reduced  pressure  of  760  torr, 

^  —  2  -6  -2  -2 

we  find  that  jjj/p0  =  9.3  •  10  (A  cm  torr  ).  According  to  Fig.  1  this  yields 

a  cathode  fall  potential  of  Vc  =  250  V,  and  pQ  dc  =  1.0  (cm  torr),  approximately. 

Therefore,  we  compute  d  =  1.3  •  10 cm,  and  E„  =  1.92  •  10^  V/cm.  The  bulk 

c  c  6  3 

energy  input  into  the  cathode  fall  region  therefore  is  j^  •  E£  =  1.04  •  10  (W/cm  ). 

Using  an  approximate  value  of  kQ  =  2,0  •  10  ^  (W/cm  K),  then  Eq.  (21)  yields 
Tc  -  TQ  =  1460  K  for  the  temperature  rise  in  the  cathode  sheath.  Assuming 
a  cooled  cathode  wall  temperature  of  about  300  K,  the  gas  temperature  in  the 
sheath  is  found  to  be  about  1760  K,  corresponding  to  a  pressure  of  about  8.8  atm 
at  constant  density.  Because  of  the  small  volume  of  the  cathode  sheath  the  char¬ 
acteristic  time  for  the  relaxation  of  discharge  generated  pressure  rise  is  much 
smaller  than  in  the  cavity.  Equilibrium  calculations  were  therefore  performed 
at  T  =  1750  K  for  several  pressures  from  8.8  atm  down  to  1  atm  initial  cavity 
pressure.  The  results  are  shown  in  Fig.  6l.  From  the  resulting  gas  composi¬ 
tion  it  is  concluded  that  enough  C02  is  dissociated  into  CO  and  O  to  consume 


Effect  of  Hydrogen  Concentration  on  Temperature  Around  Cycle 
With  and  Without  Lasing  (15th  Pulse,  PULSAR) 


p  (atm) 


Fig.  6l  -  Contaminant  Formation  in  Cathode  Sheath  (for  Equilibrium 
at  Assigned  Temperature  and  Pressure) 


almost  all  available  in  the  formation  of  HgO.  While  the  H^O  population  in 
the  cathode  fall  is  predicted  to  be  about  2  orders  of  magnitude  higher  than  in 
the  positive  column,  populations  of  OH  and  NO  are  about  the  same  as  their 
respective  peaks  reached  in  the  positive  column  at  the  end  of  the  discharge 
(see  Fig.  60).  In  contrast  to  the  main  discharge,  populations  of  NO,,  HO,  and 
N^O  do  not  even  reach  levels  of  10~  ppm.  This  is  significant  because  in  the 
absence  of  N02,  OH  and  NO  are  the  only  principal  reactants  for  the  formation 
of  HNO^.  Their  populations  being  comparable  to  those  in  the  positive  column 
it  is  concluded  that  the  cathode  fall  region  should  certainly  not  enhance  the 
formation  of  HN02.  Without  significant  amounts  of  N02>  NOj  nor  H02  which 
are  the  principal  reaction  partners  for  OH,  H20  and  NO,  respectively,  in  the 
formation  of  HNO^,  the  latter  also  has  little  chance  of  being  generated  from 
reactants  formed  in  the  cathode  fall  region  under  the  given  conditions.  Finally, 
using  Eqs.  (25)  and  (26)  it  is  found  that  when  the  H20  generated  in  the  cathode 
fall  region  is  allowed  to  mix  with  that  generated  in  die  positive  column,  the 
overall  increase  in  H20  in  the  gas  amounts  to  less  than  1%  under  the  given 
conditions. 

4.2  ARMY  CLOSED  CYCLE  GAS  CIRCULATOR 

The  Army  small  scale  closed  cycle  circulator  concept  is  discussed  in 
detail  in  Ref.  18  Its  principal  purpose  is  to  furnish  experimental  data  for 
the  design  and  development  of  high  energy  closed  cycle  laser  systems.  Some 
of  the  main  unanswered  questions  relate  to  the  closed  cycle  gas  stability, 
plasma  formation,  level  pumping  ability  and  shock  wave  attenuation  ability. 

The  problem  dealt  with  in  the  present  analysis  is  that  of  gas  stability,  par¬ 
ticularly  the  problem  of  plasma  by-product  formation  and  resulting  perfor¬ 
mance  degradation. 

Basic  operating  parameters  for  the  small  circulator  are  listed  in  Table  7. 
The  first  question  that  arises  in  setting  up  the  analysis  is  that  of  die  appropriate 
discharge  potential.  From  the  Boltzmann  code  results  for  the  COj/N-j/He  = 
1/2/3  mixture  (see  Fig.  6)  we  find  that  the  discharge  voltage  should  be  set 
such  that  an  E/N  of  about  2  *  10"  '  ^  (Vcm“)  results  which  ensures  optimum 
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Table  7 

ARMY  CLOSED  CYCLE  GAS  CIRCULATOR  OPERATING  DATA* 


Gas  Composition 

=  1/2/3 

Flow  Rate 

140  to  640  g/sec 

Cavity  Pressure 

760  torr 

Cavity  Inlet  Temperature 

200  K 

E-Beam  Energy 

125  to  200  keV 

Post-Foil  E-Beam  Current  Density 

1  to  100  mA/cm^ 

Discharge  Potential 

10  to  50  kV 

Discharge  Energy 

1000  j/lit  (max) 

Electrode  Gap 

5  cm 

Cavity  Volume 

5  x  5  x  25  cm 

Gain  Length 

25  cm 

Pulse  Width 

E -Be,.m 

2  to  30  jzsec 

Discharge 

1  to  30  jxsec 

Data  taken  from  Ref.  18. 
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upper  laser  level  pumping.  Hence,  for  given  temperature  and  pressure,  the 
desirable  electric  field  is  calculated  to  be  approximately  7.4  kV/cm.  Using 
the  given  electrode  separation  of  5  cm,  tbc  optimum  discharge  potential  be¬ 
comes  37  kV,  which  is  well  inside  the  range  indicated  in  Table  7. 

The  next  problem  which  must  be  addressed  is  the  appropriate  electron 
beam  current  density.  Assuming  an  effective  (post-foil)  electron  beam  energy 
of  100  keV,  and  pulse  widths  of  8  /xsec  and  10  jxsec  for  the  electron  beam 
the  discharge,  respectively,  the  post-foil  electron  beam  current  density  is 
limited  by  the  discharge  power  which  can  be  generated.  The  computed  depend¬ 
ence  of  discharge  energy  loading  on  electron  beam  current  density  (for  open 
cycle  operation)  is  shown  in  Fig.  62,  and  shows  that  for  a  maximum  energy 
loading  of  1  kj/lit,  the  electron  beam  current  density  should  not  exceed  40  mA/ 
cm  .  For  the  analysis  which  follows  we  have  selected  a  value  of  30  mA/cm^. 

The  application  of  our  model  to  closed  cycle  kinetics  calculations  requires 
the  gasdynamic  conditions  around  the  cycle  to  be  known  so  that  heat  exchangers 
can  be  simulated  appropriately.  Unfortunately  no  such  information  is  presently 
available  for  the  small  circulator.  As  we  had  to  recognize  previously,  in¬ 
accurate  circulator  data  can  lead  to  very  erroneous  results  (see  Fig.  50). 
Therefore,  in  order  to  perform  the  analysis  at  least  an  approximate  gas- 
dynamic  model  had  to  be  constructed  which  is  summarized  in  Table  8.  The 
table  lists  the  main  components  and  their  location  along  the  circulator  as 
interpreted  from  Fig.  4  of  Ref.  18.  Approximate  flow  cross-sectional  areas 
(denoted  A)  are  listed  also.  The  circulator  contains  three  heat  exchangers, 
the  first  two  of  which  use  water  as  coolant.  Since  the  water  is  presumably 
stored  at  room  temperature,  we  have  assumed  a  gas  exit  temperature  of  300  K 
for  these  heat  exchangers.  As  the  gas  passes  through  the  blower,  its  tempera¬ 
ture  is  assumed  to  rise  to  approximately  350  K.  Discharge  energy  input  in 
the  cavity  was  assumed  to  heat  the  gas  to  roughly  600  K,  and  the  remaining 
values  were  fortunately  known  to  be  200  K.  Knowing  the  flow  cross  sections 
and  the  temperatures,  one-dimensional  continuity  and  the  equation  of  state  can 
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Discharge  Energy  Loading  (j/lit  atm) 


Fig.  6 L  -  Discharge  Energy  Loading  as  a  Function  of  Effective 
E-Beam  Current  Density  (Open  Cycle) 


Table  8 

ARMY  CLOSED  CYCLE  GAS  CIRCULATOR  DATA 


Station 

X 

(cm) 

A 

(cm2) 

0 

Cavity  Entrance 

0 

125 

1 

Cavity  Exit 

8 

125 

2 

Diffuser  Exit 

72 

324 

3 

Duct  Exit 

160 

324 

4 

Heat  Exchanger  1 

220 

324 

5 

Blower  Exit 

290 

324 

6 

Duct  Exit 

555 

324 

7 

Heat  Exchanger  2 

630 

324 

8 

Duct  Exit 

880 

324 

9 

Heat  Exchanger  3 

975 

324 

10 

Diffuser  Exit 

1115 

390 

0 

Nozzle  Exit 

1190 

125 

u 

(m/sec) 

t 

(msec) 

31.6 

■QS 

94.6 

■M 

36.6 

1 

36.6 

18.3 

21.4 

94.0 

21.4 

220.0 

18.3 

260.0 

18.3 

390.0 

12.2 

460.0 

10.1 

580.0 

31.6 

620.0 

See  discussion  in  text. 
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be  used  to  arrive  at  an  expression  for  the  velocity  in  various  sections  of  the 
circulator,  i.e., 


u 


,m_R  JT 
1  p  '  A 


(36) 


Velocities  listed  in  Table  8  were  calculated  from  Eq.  (36)  assuming  a  mass 
flow  of  m  =  450  g/sec  and  a  constant  pressure  of  p  =  760  torr.  Finally,  die 
last  column  of  Table  8  lists  the  elapsed  time  as  computed  by  integrating  the 
reciprocal  velocity  distribution  around  the  circulator.  The  temperature  versus 
time  history  as  needed  in  the  plasma  kinetics  calculations  is  shown  in  Fig.  63. 
Note  that  the  temperature  up  to  the  entrance  of  the  first  heat  exchanger  will 
actually  be  computed  by  the  plasma  kinetics  code.  Typical  computed  (first 
cycle)  temperatures  are  shown  in  Fig.  64. 


One  of  the  questions  raised  with  regard  to  closed  cycle  operation  was 
the  effect  of  water  vapor  in  die  laser  gas  on  performance  and  contaminant 
formation.  We  have  attempted  to  answer  this  question  by  comparing  dry  gas 
calculations  with  those  assuming  10  ppm  H^O  in  the  laser  gas.  Figure  65 
shows  a  comparison  of  the  discharge  energy  loading  as  a  function  of  recircu¬ 
lation  for  the  two  cases.  Note  that  while  the  dry  gas  shows  a  slight  increase 
(about  5%)  with  recirculation  before  leveling  off,  the  moisture  contaminated 
gas  shows  a  12%  decrease  in  energy  loading  for  the  first  pulse,  which  then 
decreases  further  with  recirculation.  Corresponding  first  pulse  discharge 
current  density  distributions  are  shown  in  Fig.  66.  When  averaged  over  the 
discharge  pulse  width  of  10  jusec  these  current  density  distributions  represent 
average  discharge  currents  of  1478  and  1304  A,  respectively. 

It  is  surprising  to  see  how  little  effect  the  differences  in  discharge 
energy  loading  have  on  the  small  signal  gain  distributions  which  are  shown 
in  Fig.  67.  While  the  peak  gain  remains  virtually  the  same  (at  12%  cm  ') 
the  only  effect  is  a  slight  shift  of  the  peak  with  time.  It  should  be  remembered, 
however,  that  the  analysis  of  the  PULSAR  showed  the  same  phenomenon.  While 
for  that  device  small  signal  gain  distributions  were  almost  identical  (see  Fig.  46), 
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Fig.  63  -  Temperature  Distribution  Around  Circulator  as  a  Function  of  Time 


Fig.  65  -  Discharge  Energy  Loading  as  a  Function 
of  Recirculation  and  H^O  Contamination 
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C°2/N2/He  =  1/2/3  (CCGC) 
E/N  =  2.0  •  10' 16  Vcm2 
j£B  =  30  mA/cm2 
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Fig.  66  •  Effect  of  Moisture  on  First  Pulse  Discharge 
Current  Density 
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the  corresponding  pulse  output  energies  differed  substantially  (see  Fig.  55).  It 
must  be  concluded  that  small  signal  gain  distributions  by  themselves  are  not 
necessarily  good  indicators  for  the  performance  of  a  device. 


The  decrease  in  discharge  energy  loading  due  to  moisture  in  the  gas  can 
be  explained  by  examining  by-product  concentrations  in  the  gas  as  a  function 
of  recirculation.  In  a  gas  with  no  trace  of  hydrogen,  the  various  oxides  of 
nitrogen  really  have  no  major  effect  as  long  as  only  small  quantities  are 


created.  Figure  68  shows  concentrations  of  CO,  and  O^,  as  well  as  of  NO, 
NO^  and  N^O  as  they  develop  during  recirculation  of  a  perfectly  clean  CO^/N^ 
He  =  1/2/3  mixture.  These  concentrations  should  be  compared  to  those  shown 


in  Fig.  69  which  result  when  the  gas  is  moisture  contaminated.  The  most 


noticeable  difference  is  that,  for  the  moist  gas,  the  resulting  NO^  concentra¬ 
tion  is  roughly  an  order  of  magnitude  lower  than  that  for  the  dry  gas.  Except 


for  03,  which  for  the  moist  gas  shows  not  only  a  much  lower  concentration  but 
also  a  reversal  in  trend,  the  other  species  do  not  show  any  drastic  change. 


The  principal  reason  for  die  decrease  of  the  NO^  concentration  is  that 
N02  reacts  with  H20  to  form  HN02  and  OH.  N02  further  reacts  with  OH  to 
form  HNO^  while  NO  combines  with  the  remaining  OH  radicals  to  generate 
more  HN02<  Figure  70  shows  that  roughly  90%  of  the  available  H20  is  dis¬ 
sociated  to  formH^.  The  remaining  10%  (1  ppm  in  this  case),  however,  is 
sufficient  to  combine  with  the  available  1  ppm  NO^  (in  the  dry  gas)  to  start 
the  chain  reaction  of  forming  HN02  and  HNO^.  Under  these  conditions, 
ignoring  the  NO  for  the  moment,  we  would  expect  to  see  roughly  1  ppm  of 
HN02/HN03  combined  being  created.  Figure  70  shows  that  this  is  precisely 
the  case.  We  also  conclude  that  formation  of  HO^NO^  plays  only  a  very 
minor  role.  It  is  generated  mainly  from  HO,  and  N,Oc  both  of  which  are 

u  CD 

present  only  in  negligible  amounts  under  the  present  conditions. 


'  li  CO^/N^/Hr  1/2/3  (CCGC) 

¥  E/N  =  2.0  •  10'16  Vcm* 

I  jEB  =  30 
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Fig.  68  -  Dry  Gas  Contaminant  Generation  During  Closed  Cycle  Operation 
(Without  Lasing) 


Cycle  Number 


Fig.  69  -  Moist  Gas  Contaminant  Generation  During  Closed 
Cycle  Operation  (Without  Lasing) 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  overall  objective  of  the  effort  described  in  this  report,  namely,  to 
extend  the  existing  EDLAMP  computer  program  to  account  for  the  electron 
energy  distribution  in  the  plasma  chemistry  analysis  of  closed  cycle  EDL 
systems,  and  to  evaluate  specific  operational  gas  mixtures  has  been  accom¬ 
plished.  The  coupling  of  a  Boltzmann  code  into  the  analysis  via  generalized 
curve  fits  of  the  drift  field  velocity  and  the  various  electron  impact  excitation 
rates  has  added  substantially  to  the  usefulness  of  the  code  in  analyzing  CO^ 
electric  discharge  lasers  using  gas  mixtures  consisting  of  CO^*  N^,  He,  CO 
and  in  arbitrary  proportions.  Implementation  of  several  other  refinements 
such  as  a  revised  model  of  electron-beam  ionization,  preparation  of  a  set  of 
updated  general  reaction  mechanisms  and  rates,  and  incorporation  of  at  least 
a  simple  model  to  assess  the  effect  of  cathode  sheath  conditions  on  the  plasma 
chemistry  in  the  discharge  have  further  improved  the  model  as  an  analytical 
tool.  The  improved  model  has  been  tested  by  comparing  our  predictions  with 
supplied  experimental  results  in  one  case,  and  by  comparison  of  our  analyt¬ 
ical  results  with  those  of  others  in  another  case.  The  good  agreement  achieved 
lends  credibility  to  our  theoretical  model  and  its  application  to  the  analysis  of 
closed  cycle  EDL  systems. 

A  detailed  analysis  has  been  presented  for  two  closed  cycle  EDL  de¬ 
vices  which  are  presently  under  development.  For  the  PULSAR  device,  which 
uses  a  CO^/N^/H^  =  1/3/0.08  gas  mixture,  the  effects  of  variations  in  dis¬ 
charge  voltage,  electron-beam  current  and  concentration  on  laser  pulse 
energy  has  been  investigated.  Variations  in  the  first  two  parameters  basic¬ 
ally  served  to  define  optimum  operating  conditions  subject  to  hardware  con¬ 
straints  (maximum  discharge  energy  loading).  Variations  inH^  content  were 
studied  with  the  aim  of  minimizing  gas  contamination  under  closed  cycle  oper¬ 
ation.  It  was  found  that  decreasing  the  concentration  actually  increases 
the  concentrations  of  electrophilic  species  formed  (such  as  HNO^  and  HNO^) 
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which  results  in  decreased  output  pulse  energy  as  well  as  a  stronger  rate 
of  decay  under  closed  cycle  operation.  From  an  assessment  of  the  conditions 
in  the  cathode  sheath  for  this  device  it  is  concluded  at  this  time  that  the 
cathode  fall  region  should  not  significantly  contribute  to  the  contamination 
problem. 

Analysis  of  the  Army  circulator  device  centered  around  studying  the 
effect  of  contamination  of  the  CO^/N^/He  =  1/2/3  gas  mixture  by  water 
vapor  on  closed  cycle  discharge  energy  loading  and  small  signal  gain.  While 
the  effect  on  small  signal  gain  was  found  to  be  negligible,  the  discharge  energy 
loading  decreased  noticeably  due  to  the  formation  of  HNO^  and  HNO 

A  series  of  calculations  was  also  performed  to  study  the  effect  of  con¬ 
taminant  removal  via  catalytic  reactor.  This  was  done  for  the  CO^/CO/N^/He  = 
1/0.25/5.5/0.125  gas  mixture  in  the  Army  closed  cycle  circulator,  assuming 
a  very  simplified  temperature  history  around  the  cycle.  These  calculations 
were  completed  before  the  importance  of  an  accurate  temperature  history  was 
recognized  during  the  analysis  of  the  PULSAR  device.  Since  the  results  must 
therefore  now  be  considered  in  error  these  calculations  have  not  been  docu¬ 
mented  here.  Unfortunately  time  did  not  permit  repeat  of  this  analysis  during 
the  present  effort.  A  preliminary  conclusion  reached  from  analyzing  various 
discharge  conditions  for  a  variety  of  gas  mixtures  is  that  the  selection  of  dis¬ 
charge  voltage  and  electron-beam  current  density  such  that  optimum  efficiency 
is  achieved  for  the  given  maximum  discharge  energy  loading  may  be  one  of  the 
best  ways  to  minimize  the  closed  cycle  discharge  generated  contaminant  problem. 

The  comparison  of  the  code  results  with  experimental  data  and  its  appli¬ 
cation  to  analyzing  test  cases  has  brought  several  specific  deficiencies  to  light 
which  should  be  remedied  in  a  future  effort. 

1.  While  we  do  not  know  the  reason  for  the  slight  discrepancies  be¬ 
tween  our  analytical  and  Freed's  experimental  results,  it  may  be 
possible  to  improve  the  agreement  by  further  refinement  of  the 
kinetics  model.  As  new  kinetics  information  appears  in  the  litera¬ 
ture,  it  must  be  extracted,  evaluated  and  implemented.  This  re¬ 
quires  an  almost  continual  effort  to  keep  abreast  of  the  latest 
developments  in  kinetics  research. 
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2.  Whereas  our  model  computes  small  signal  gain  values  in  good 
agreement  with  experiment  for  low  discharge  energy  loading 
(Freed's  case,  about  0.1  kj/lit  atm),  predicted  small  signal 
gain  for  high  energy  loading  (ABEL,  about  0.8  kj/lit  atm)  are 
judged  too  high  although  the  output  power  appears  to  be  approxi¬ 
mately  correct.  A  possible  reason  for  this  behavior  of  the  model 
may  be  the  truncation  of  COg  vibrational  levels  above  the  CC>2(001) 
level.  This  simplification  may  not  be  valid  for  high  energy  loading 
situations  accompanied  by  comparatively  higher  gas  temperatures 
leading  to  non-negligible  populations  of  the  ignored  higher  energy 
levels.  Further  analysis  of  high  energy  loading  cases  and  com¬ 
parison  with  reliable  small  signal  gain  measurements  should  be 
undertaken  to  resolve  this  question. 

3.  Regarding  the  EDLAMP  code  and  its  capability  to  predict  closed 
cycle  performance  the  most  pressing  need  is  a  revision  of  the 
catalytic  reactor  model  option.  Provision  should  be  made  to 
"scrub"  certain  species  not  at  the  end  of  the  cycle  (as  presently 
done)  but  somewhere  between  cavity  exit  and  the  first  heat  ex¬ 
changer  because  of  the  high  temperature  requirements  for  effi¬ 
cient  catalyst  operation.  Furthermore,  contaminant  removal 
efficiency  should  be  redefined  to  reflect  actual  catalyst  capability 
to  remove  certain  species  down  to,  for  example,  one  part  per 
million  concentration. 

4.  An  effort  should  be  undertaken  to  incorporate  a  cathode  sheath 
model  in  the  EDLAMP  code  in  a  more  rigorous  fashion  than 
what  was  possible  under  the  present  effort. 

5.  Further  analysis  of  test  cases  should  include  the  effect  of  varia¬ 
tions  in  resonator  output  coupling  on  closed  cycle  laser  perform¬ 
ance  (results  presented  here  have  arbitrarily  assumed  an  output 
coupling  of  50%),  as  well  as  other  combinations  of  discharge  voltage 
and  electron-beam  current  density  (for  given  maximum  discharge 
energy  loading). 

6.  Closed  cycle  laser  output  performance  predictions  should  be  com¬ 
pared  with  experimental  results  as  these  become  available.  Of 
particular  interest  is  the  validation  of  predicted  decay  rates,  with 
and  without  catalytic  removal  of  key  contaminant  species. 

The  iast  recommendation  made  is,  of  course,  intimately  related  to 
electrophilic  plasma  byproduct  generation.  One  of  the  key  experiments  that 
can  thus  be  suggested  should  be  the  measurement  of  NO,  N^O,  HNO.,,  HNO^ 
and  H^O  concentrations.  The  analysis  has  shown  that  these  measurements 
should  be  taken  somewhere  just  upstream  of  the  cavity  entrance.  Measure¬ 
ments  of  OH  and  NO^  concentrations,  if  possible,  should  be  attempted  in  die 
cavity  just  after  the  discharge  pulse  has  terminated. 
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Appendix 

REACTION  RATE  COEFFICIENTS 


A  Ionization  Reactions 

B  Attachment  Reactions 

C  Detachment  Reactions 

D  Ion-Electron  Recombination  Reactions 

E  Ion -Ion  Recombination  Reactions 

F  Charge  Exchange  Reactions 

G  Electron  Impact  Excitation  Reactions 

H  Vibrational  Relaxation  Reactions 

J  Ternary  Recombination  Reactions 

K  Binary  Reactions 


NOTE:  Reactions  A1  through  A19  are  of  type  14  or  15  (see  Vol.II,  EDLAMP  User's  Manual). 
For  program  uo”,  only  the  constant  factor  of  the  rate  coefficient  is  read  (the  factor 
-0  23 

£EB  coded  into  the  program). 


I  £ 


O  •— «  (Nl  •-* 


OcOOO^OO^m 


4> 

4) 

a 

e 

u 

e 

a) 

d 

d 

N 

d 

\ 

4) 

d 

N. 

\ 

in 

m 

\ 

O' 

• 

fM 

• 

o 

• 

o 

• 

o 

• 

in 

• 

<M 

m 

co 

CO 

CO 

CO 

H 

H 

i 

i 

i 

H 

| 

| 

* 

I 

CL 

& 

9- 

a 

a 

X 

X 

X 

a 

X 

X 

a> 

0) 

V 

X 

V 

4) 

m 

in 

4) 

00  00 
I  I 


o  o 

•  • 

m  in 


•  O  vO  i 

(NJ  *  * 

i  <u  <M  N  j 

n  I  4)  I  4> 
X.  d  a 
o  ^  ^ 


*— 4  O' 

I  I 


o  o 
•  • 
in  in 


•-«  fM 

H 

I  I 

o  o 


•  l. 


/•'-N  ' 


o 

X 

X 

p-4 

pH 

m 

O 

X 

X 

o 

o 

o 

2 

2 

2 

O 

+ 

+ 

+ 

+ 

z 

z 

+ 

+ 

+ 

+ 

1 

•  (M 

i  eg 

I  ro 

+ 

1  (Vi 

+ 

1 

1  fM 

1  <M 

■ 

X 

o 

o 

o 

1 

o 

1 

O 

o 

o 

o 

o 

z 

z 

z 

o 

z 

o 

1  I 


H  N 

CQ  «  CQ 


z  z 
+  + 


—  (NJ 

E  O 
<M  Z 


eg 

O  O  tM 


I  CO  + 

o  . 


I  I  t  I  t  t  I  I  1  I  I  I 


H 

m 

<1) 

H 

% 

s 

+ 

4) 

2 

4) 

<U 

4) 

+ 

V 

+ 

V 

+ 

O 

O 

o 

a> 

+ 

+ 

CM 

+ 

co 

+ 

'j* 

V 

+  e 

+  e 

V 

+ 

+ 

in 

4) 

+ 

fM 

H 

+ 

+ 

fM 

o 

o 

0 

0 

+ 

(M 

o 

o 

+ 

+ 

X 

+ 

<M 

eg 

O 

N 

2 

z 

z 

0 

0 

0 

eg 

<M 

(M 

in 

o 

o 

0 

u 

53 

X 

53 

53 

z 

z 

z 

Z 

z 

0 

O 

U 

O  f-H  CM  cO  in  sD 

^  in  vO  r-  00  O'  H  p4  H  H  H  H 

ffl  cqpQcQCQ«mcQCQcQCQCQPQ 


.  i 


m  m 

sT  M 
•  •  • 

D  O  O 

l  I  i 

H  H  H  _  _  ^  _  _  _  _ 

oo 

0s  0s  00  H  pH  0s  H  0s  H  H  pH 

I  I  I  I  I  I  I  I  I  I  I  I 


oomoocNJO^oooin 


co  o  o 

•  •  • 

cO  ^ 


<U  <U 

+  + 


u 

o 

fM 

O 

V 

+ 

4) 

+ 

V 

<v 

V 

+ 

i) 

+ 

<u 

+ 

+  e 

<u 

fM 

o 

co 

O 

fM 

o 

(NJ 

fM 

+ 

+ 

O 

fM 

fM 

fM 

fM 

+ 

+ 

+ 

o 

o 

O 

fM 

o 

o 

O 

cO 

fM 

co 

j  1 

o 

a 

2 

z 

o 

X 

X 

EC 

Z 

o 

o 

o 

I  I  I  t  I  I  I  t  1  I  II 


o 

O 

o 

fM 

EC 

0 

z 

z 

o 

H  + 

o  + 

EC 

Z 

0 

o 

CO 

o 

+ 

+ 

+ 

+ 

+ 

( 

1 

+ 

+ 

+ 

+ 

+ 

1 

i 

1 

i 

1 

EC 

EC 

1  fM 

1  fM 

I  fM 

1  CO 

1  CO 

o 

0 

0 

0 

0 

0 

o 

o 

0 

0 

0 

o 

^H 

fM 

on 

in 

o 

r- 

00 

O' 

o 

H 

pH 

pH 

fM 

H 

U 

o 

o 

u 

O 

U 

U 

U 

u 

U 

U 

U 

\P 


•  r 


VS*.  O'O'O'O'O'C'O'O'O'O'O'O'O'O'O'O' 
“nnnnnnnnnnnnnnnn 

as 


u  r- 
C3  ‘ 


o  o 

m  m 


o  o  o  o  o 

•  •  •  I  • 

m  co  co  co 


(VJ 

o  o 


o  o 

(Vj  (Vj 


oo 

O 

(VJ 

O 

u 

+ 

u 

+ 

O 

+ 

CVj 

o 

0 

a 

(M 

o 

oo 

0 

+ 

(VJ 

o 

(VJ 

o 

Ps] 

U 

O 

oo 

po 

+ 

o 

o 

u 

+ 

o 

U 

o 

53 

X 

+ 

+ 

O 

0 

PO 

4- 

4- 

4- 

4- 

+ 

4- 

+ 

+ 

Z 

z 

0 

—v 

o 

o 

PO 

o 

(vj 

+ 

+ 

+ 

o 

(VJ 

o 

(Vj 

o 

(VJ 

+ 

o 

(VI 

o 

(VI 

O 

(VJ 

O 

(VJ 

o 

00 

O 

O 

<— ** 

O 

O 

o 

O 

O 

o 

o 

o 

H 

H 

I— < 

' 

" 

fH 

PO 

(v| 

(vj 

(VJ 

(VJ 

'w” 

(M 

00 

(VJ 

'Ini 

O 

O 

O 

0 

0 

O 

0 

o 

0 

O 

O 

o 

O 

O 

U 

O 

u 

U 

u 

o 

o 

o 

u 

u 

U 

o 

o 

U 

nO  vO  nO  vO  nD 

I  I  I  I  I 


rsl  TT  Tf  rr  ^ 

•  ••••• 
H  CVj  (Vl  <NJ  (Vj  <\| 


o  o  + 

PO  (M  (M 

33  33  O 


t  t  It  t  It  I  t  t  t  I  ft  It  1 


1 

i  co 

I  TT 

1  (VJ 

»  co 

1  CO 
0 

i 

0 

1  P0 

i  po 
0 

t  Tf 
0 

I  00 

1  PO 

o 

i  'i* 

o 

i  (VJ 

o 

1  PO 

o 

1  (VJ 

0 

0 

0 

o 

I  P0 

I  -0< 

u 

U 

0 

U 

u 

0 

o 

u 

z 

2 

o 

L 

u 

u 

z 

z 

0 

0 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

4- 

+ 

+ 

r 

» 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+  04 

+  P0 

+  PO 

+ 

+ 

4* 

+ 

+ 

4-  (vj 

+  00 

+  PO 

+  P0 

+  PO 

+  PO 

0 

0 

0 

0 

0 

O 

o 

o 

o 

o 

o 

0 

0 

0 

0 

0 

0 

u 

U 

u 

u 

u 

o 

po 

PO 

CO 

PO 

PO 

l 

u 

u 

u 

U 

u 

O 

33 

33 

33 

33 

X 

X 

X 

X 

5; 

X 

X 

o 

•— < 

00 

PO 

'T 

m 

vO 

r- 

1 

H 

PO 

m 

in 

o 

r- 

00 

O' 

H 

r—« 

1—4 

H 

i-4 

w 

w 

W 

w 

W 

W 

W 

W 

W 

W 

W 

w 

W 

w 

W 

W 

W 

CJNOOOOCT'O'O'OOOCT' 

fifsjtvirOrOtOfslMMroc^fOrJ 

« 


o  o  o  o 

•  •  •  • 

sO  so  «n  CO 


o  o  o 

*  •  • 

r<">  oo 


4 

4 

4* 

4 

0 

O 

(M 

CM 

z 

£ 

O 

0 

CM 

o 

o 

O 

o 

0 

4 

4 

4* 

4 

o 

4 

4* 

4 

+ 

4 

o 

o 

fM 

O 

o 

fM 

o 

£ 

+ 

z 

4 

4 

z 

fM 

o 

o 

(M 

O 

O 

(M 

O 

o 

4 

o 

4 

4 

o 

(M 

fM 

fM 

fO 

4 

4 

fM 

fM 

(M 

<f> 

4 

O 

O 

o 

0 

(M 

o 

o 

O 

o 

o 

(V 

u 

u 

z 

z 

o 

z 

o 

o 

z 

z 

o 

!  t  t  t  t  l  t  1  i  i  t  f  l 


i  m 

o 

1  M* 

o 

1  fM 

o 

i  co 
0 

*  <M 

1  M* 

t  41 

i  M* 

1  fM 

1  (O 

u 

o 

z 

z 

o 

0 

0 

0 

o 

0 

i  fM 

4 

+ 

4 

4 

4 

4 

o 

u 

z 

z 

o 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

0 

0 

0 

0 

0 

0 

4  fM 

4  fM 

4  fM 

4-  fM 

4  fM 

z 

z 

z 

z 

z 

z 

0 

o 

0 

0 

0 

o 

H 

fM 

CO 

in 

vO 

00 

O' 

fM 

CM 

(M 

fM 

fM 

CM 

CM 

fM 

fM 

fM 

w 

w 

w 

w 

w 

w 

w 

w 

w 

« 

!  ii 

;  i; 


gj  co  co  co  co  co  ^ 


COCO<OCOCO»-HCOCOCOCO 


* 

! 

V 

, _ _ 

_ _ _ 

, — „ 

J 

4-1 

O 

o 

. — „ 

^-v 

o 

o 

. — . 

, — _ 

. - s 

o 

o 

o 

« 

•  . 

H 

■-H 

O' 

H 

H 

O' 

O' 

r-H 

o 

O' 

O' 

o' 

O 

pH 

pH 

pH 

oj 

1 

1 

1 

1 

1 

1 

» 

1 

1 

1 

1 

1 

*H 

I 

1 

1 

• 

x— ' 

x»p 

x— -- 

x«w» 

w 

X—* 

”x»p 

X—" 

X— * 

•x-p 

X— ' 

•x.P’ 

x^- 

x— p 

x—o» 

*  * 

v£> 

o 

M< 

TT 

LO 

r-H 

OJ 

O 

rsj 

O 

o 

o 

04 

CO 

IT) 

m 

NO 

oi 

xO 

H 

H 

TT 

r4 

r-H 

co 

pH 

CO 

co 

cO 

pH 

04 

'j* 

04 

O 

o 

o 

iaj 

o 

<u 

o 

^OJ 

o 

o 

o 

r*~4 

o 

o 

~1\] 

o 

i—4 

o 

^OJ 

o 

E 

O 

o 

X 

o 

o 

(M 

0 

u 

X 

u 

u 

u 

4- 

0 

o 

o 

+ 

4- 

4- 

u 

0 

o 

U 

4- 

u 

4- 

+ 

4- 

O 

z 

z 

z 

4- 

4-  <M 

4-  (M 

+ 

4- 

4- 

4- 

4- 

4- 

4- 

4- 

4- 

+ 

4- 

+ 

o 

o 

o 

+ 

4- 

4-  CM 

o 

4- 

o 

o 

o 

4- 

4- 

4-  (M 

4- 

+  CM 

u 

o 

u 

o 

+  <M 

0 

CO 

o 

u 

CO 

CO 

o 

o 

o 

o 

O 

E 

X 

X 

2 

o 

u 

E 

z 

X 

X 

E 

o 

u 

o 

u 

It  It  I  I  It  t  I  II  t  tt  M 


o 

o 

o 

o 

o 

o 

o 

o 

0 

w ’ 

o 

o 

o 

o 

o 

"cm 

(VJ 

0 

o 

cm 

(M 

OJ 

o 

o 

(M 

CM 

O 

CM 

0 

E 

Z 

u 

E 

E 

o 

^04 

CM 

o 

O 

E 

E 

E 

z 

O 

U 

4- 

+ 

+ 

4- 

4- 

U 

0 

o 

o 

CM 

+ 

4- 

4- 

4- 

+ 

4- 

+ 

4- 

+ 

4-  CM 

4-  (M 

4-  M 

+ 

u 

O 

u 

O 

4-  CM 

4-  CM 

4-  CM 

4-  CM 

+  CM 

4-  cm 

4- 

0 

0 

O 

0 

0 

+ 

+ 

4- 

4- 

4- 

O 

0 

0 

0 

0 

0 

o 

u 

u 

o 

O 

V 

V 

4- 

4- 

4- 

4- 

u 

u 

o 

u 

<J 

U 

o 

E 

E 

E 

E 

E 

X 

z 

z 

z 

z 

o 

pH 

(M 

CO 

T 

in 

sO 

r~ 

»"H 

(M 

CO 

rr 

m 

vO 

(X 

oo 

O' 

rH 

pH 

pH 

l—H 

pH 

1*4 

E 

E 

E 

E 

E 

E 

E 

E 

U i 

1*4 

t*4 

1*4 

\*4 

u« 

1*4 

co  ro  •-<  co  co  co 


00  O'  O'  O'  O'*  O' 

co  co  co  co  co  co 


O  >-s  O  CM 

H  O  H  H 

I  I  I  I 


OOOO  0^000 

•— «  «— •  »-4  •-*  H  (JV  H  H  H 

I  I  I  I  I  I  I  I  I 


Xr^^o^cMO^oom 

O  CO  CO  0s  H  H  H  ^ 


OOOO 
•  •  *  • 

CO  co  co  CO 


H  CO  N 


03  w' 

2  0  2  O 

2  +  U  +  O  + 

+  +  <M  +  +  +  + 

+  <M  O  4-  (M  O  +  (M  O 
0  0  0  2  0  2 


1  CO  I  CO 

o  o 

2  O 


I  tt  I  t  t  I  I  I  I  f  t  ft  t  t 


(M  (M  O 

o  o  o 


<M  (M 

)  o 


o  o 

+  + 

+•  N  + 

2  O 


(VI  <VJ 

2  O 


+  +  1-  -r  4. 

+  +  +  <v)  #  <M  +  +  # 

+  +<M+(vlO  O  *N»N  O 

000UU20U 


CO 

o 

2 

E  2 


i  co  i  fo  i  vr 

O  O  O 

U  O  O 


00 

O' 

© 

H 

CM 

CO 

in 

t". 

00 

O' 

o 

H 

(M 

co 

H 

H 

CM 

CM 

03 

03 

CM 

CM 

CM 

(VJ 

CM 

CM 

CO 

CO 

CO 

CO 

to 

u* 

Ui 

( M 

Ui 

Ui 

Ui 

Ui 

Ui 

to 

to 

to 

to 

to 

to 

Ui 

v  • 


-  10 

'  !  i. 


[s.  (Nj  f''-  H  M  N  00  N 


•« 4 

u 
•— « 
v» 

S-l 

4) 

O 

O 

1 

H 

in 

CM 

i 

i 

H 

oo 

<M 

1 

1 

H 

QQ 

cm 

i 

o 

CO 

1 

O' 

CM 

1 

i 

H 

00 

CM 

l 

1 

H 

00 

(M 

1 

1 

H 

P 

CM 

1 

H 

1 

o 

H 

1 

H 

H 

1 

H 

H 

1 

o 

pH 

1 

o 

H 

1 

o 

H 

1 

V 

o 

o 

o 

o 

CM 

O 

O 

o 

00 

o 

o 

• 

00 

• 

m 

• 

CO 

• 

r- 

• 

<« 

m 

• 

co 

• 

cO 

• 

CM 

in 

CO 

CO 

CM 

pH 

CM 

00 

M4 

H 

H 

CO 

^  O 

(Jv  H 

I  I 

vO  o 


I  M  l  N 

o  o 


in 

H 

I'M 

o 

1  M 

o 

""oj 

o 

CM 

pH 

in 

2 

2 

H 

r— 4 

u 

u 

0 

u 

O 

2 

2 

+ 

+ 

2 

2 

4- 

+ 

o 

H- 

+ 

+ 

i  -r 

l  co 

+ 

+ 

1  CM 

I  CO 

~b 

i  co 

1  CO 

1  CO 

1  CO 

o 

0 

1  cO 

i  m* 

o 

o 

1  CM 

o 

o 

o 

0 

u 

2 

o 

o 

2 

2 

O 

2 

o 

+  + 

1  co  i  co 

o  o 

2  2 


ft  I  It  II  I  t  It  t  t  It  1  t 


in 

2 

o 

o 

o 

pH 

2 

pH 

pH 

2 

CM 

+ 

'~CM 

+ 

+ 

0 

o 

CM 

0 

0 

+ 

CM 

2 

2 

0 

U 

2 

0 

0 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

1  CO 

1  CO 

1 

l  CM 

1  CM 

•  CM 

1  CM 

0 

0 

0 

0 

O 

0 

0 

U 

u 

+ 

+ 

+ 

1  CO 

I  M* 

I  M* 

0 

o 

o 

u 

u 

o 

CO 

O  E 


l  M*  I  N  I  M  <  CM 

o  o  o  o 

U  2  2  2 


00 

O' 

o 

pH 

CM 

CO 

V 

in 

sO 

co 

CO 

CO 

M* 

M4 

M< 

Vf 

M4 

M1 

M* 

1*4 

Ui 

1*4 

u> 

Uc 

tn 

U4 

Uc 

Ui 

Uc 

Ui 

/  v 


}• 


d 


!  f 


it  ( 


I  I 


u 


4> 

U 

H 

o 

„ — „ 

. — * 

o 

o 

o 

o 

o 

o 

o 

o 

•H 

r-4 

O' 

o 

o 

o 

O' 

O' 

H 

H 

H 

H 

H 

H 

H 

H 

m 

1 

1 

1 

1 

| 

i 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

4) 

» 

w 

w 

'w» 

*»■** 

•w* 

O 

00 

M4 

o 

CM 

CO 

00 

00 

o 

o 

o 

M4 

o 

00 

o 

00 

in 

U 

<d 

4-> 

<4 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

•H 

H 

CO 

H 

m 

r— 4 

00 

rH 

vO 

co 

nO 

00 

H 

00 

H 

M4 

00 

co  33  33 

o  o  o 


33 

O 

+ 


CO 

O  CO  CO  CO  CO 

33  O  O  co  O  O  O 

+  4-  +  O  +  4-  + 


co  co  ro 

O  co  O  O 

+  0  +  4. 


i  Cl  l  I  Cl  I  N  +  l  CO  •  CO  I  CO  I  **■  I  Cl  I  I  I  rt  +  I  T  I  Cl 

O  33  O  O  i  co  O  O  O  33  i  co  O  33  O  O  ■  co  O  O 

zozzozzzoooozzooz 


a 

o 


t  t  I  t  I  t  !  It  t  I  M  I  t  M 


o 

rt 

o> 

oi 


o 

o 


o 

o 


CO 

o 

CO 

o 

co 

CO 

0 

CO 

0 

co 

o 

co 

co 

o 

CO 

o 

00 

O 

CO 

Z 

z 

Q 

33 

CO 

O 

o 

O 

o 

0 

+ 

E 

§ 

z 

0 

+ 

33 

Z 

0 

0 

U 

33 

z 

Z 

0 

O 

z 

i  co 

+ 

+ 

+ 

+ 

■ 

+ 

+ 

+ 

+ 

+ 

4* 

+ 

+ 

+ 

+ 

+ 

0 

i 

1 

1 

1 

33 

I  co 

l  CO 

i  co 

1  CO 

1  Cl 

*  CO 

1  Cl 

1  Cl 

1  Cl 

i  ^ 

1  V 

z 

0 

o 

0 

o 

0 

O 

0 

0 

O 

0 

0 

o 

o 

o 

o 

o 

H 

00 

CO 

'T 

in 

iO 

r~ 

oo 

O' 

o 

H 

CO 

Cl 

'f 

in 

nO 

f- 

in 

in 

in 

in 

m 

in 

in 

in 

in 

vO 

nO 

03 

vO 

vO 

iO 

vO 

o 

Ui 

U* 

U A 

Ph 

Ui 

Ui 

Ui 

Pm 

u< 

U* 

k 

Pm 

p4 

(*4 

u« 

P^ 

A- 10 


j 


H  0s  0s  0s  O' 


fl) 


pH 

s 

in 


o 

+ 

<VJ 

(M 

N 

H 

H 

H 

X 

O 

£ 

X 

in 

(Nl 

in 

tn 

U 

o 

2 

U 

O 

+ 

+ 

+ 

+ 

+ 

1  (Nj 

i  cO 

■ 

i 

1  <M 

0 

o 

o 

O 

o 

X 

1 

t 

o 

t 

t 

t 

(M 

<M 

OJ 

0 

z 

o 

o 

z 

+ 

+ 

+ 

+ 

o 

(Nl 

CM 

(Nl 

1  —t 

1  H 

1  H 

1  r-4 

+ 

X 

X 

X 

X 

1 

If) 

in 

m 

in 

0 

o 

u 

U 

U 

00 

O' 

o 

H 

(M 

o 

NO 

t- 

Is- 

r- 

fa 

fa 

fa 

Ui 

! 

A-ll 


/  '*  • 


CO  <*> 
(M  CO 


vD 

CO  CO 
•  • 

o  o 

V  V 

a  a 


j,  =  « 

3  »o  a 
J  \ 
O  O  id 
•  r*  cj 

SM* 
B1  *  a 

V  IT)  x 
^  1  « 
M  flj  « 
—l  .O 
I  H  I 


*  Q, 

2  * 

rL  V 


v  _ 

a  v 

>  a 

a  > 

2  2  _ 

CL  *  ft 

x  a  3 

o  «  K 

«  3  ^ 

a  «n  — 

r-  co  g* 

f"  «  X 

i  in  o 


•*  +  v  + 

m  Uh  co  — « 


-i  1+  o 


o  o  o 

-H  O  CO 

O  — I  o 

^o  'co  'co 

O  O  O 

O  O  O 


o  cm 

CO  o 

o  o 

'co  'co 

O  O 

O  U 


1  t!  u  II  ll  11 


co  co  co 
0  0  0 

o  o  u 


cl  co  CO 

O  O  O 


CO  co 

o  o 

O  U 


ft) 

+ 

o 

o 

a) 

+ 

0) 

+ 

H 

s 

o 

o 

<M 

o 

o 

+ 

*  CO 

0 

O 

'co 

o 

0 

U 

U 

z 

O 

U 

O  CM 

o  o 


A  - 12 


i 

'l 

l  F: 

|i 

In 

i 

i  , 

1 1. 


*Jinin»n«niniomin 

K 


O  —I  ^  O  —t  t-( 


O  ~ 

u  2, 


O  o 
O  O 


o  S  s  +  u  o 

+  +  +  o  +  O 


5  i 


O  O  0'n  +  n  O 

ooozzzzz 


1  t  It  t  I  t  tt 


<M  N 

o  o 


u  S  - 


<M  PJ  (SI 

D  O  O 


M  w  V  V 

+  2  o  u  u 


*  <M  *  *  +  +  +  +  + 

O  O  0#<'j*N*M#cg*rj 
OOOZZZZZ 


cgco^'inof'OOCT' 

H  •— 4  «— 4  H  i—4  *— 4  H  —4 

oooooooo 


A- 13 


I  j. 

s  i: 


00  00  O' 
in  m  in 


^  rsi  m  m 
4}  >0  vO  vO 


00  00 

>0  N0 


^  to 

to  to  to 

m  •  m 

•  o  • 


H  o 

o<  H 


^  2 


cm  cm 

H  H 
I  I 

?r  y 


o  o 
H  H 


i  i 


o 

in 

w  in  >n 

0-  *  * 

kT  «-l  <M 

5  H  H 

*H  00  pH 

pH  pH  CM 

I  I  I 


I  (M  pH  pH  CO 

H  H  H  H  H 

o  ^  t^-  In  S'  ^ 

H  (M  H  r-l  ^  {VJ 

•  •  I  I  I  I 


•  o  o 

PI  •  • 

H  en  p-« 

«  H  ' 

H 

im  in  /> 

l  P4  O' 

I  I 

in  %—  ■ _ ■ 

<M  N  'P 


f 

4m 

o 

• 

00 

• 

sO 

H 

• 

H 

ro 

CM 

pH 

• 

pH 

• 

m 

• 

H 

ni 

• 

■<r 

• 

• 

pH 

• 

vO 

• 

H 

h-i 

4) 

T 

O 

^ 

p-** 

u 

o 

o 

o 

£ 

pH 

H 

o 

p- 

r 

4) 

£ 

H 

H 

CM 

pH 

'**p* 

0 

o 

"(VI 

O 

o 

(M 

o 

c 

(O 

o 

CM 

CM 

m 

<M 

0 

CM 

o 

^ni 

O 

•<r 

•v 

-r 

o 

O 

l 

2 

2 

U 

O 

U 

H 

2 

2 

2 

X 

O 

2 

2 

2 

O 

mm 

4- 

+ 

+ 

*4- 

+ 

O 

O 

+ 

+ 

+ 

■4 

+ 

4- 

4- 

4- 

4- 

4- 

, 

p«* 

p.* 

p-». 

O 

O 

p-N. 

p-*» 

P-N 

«p*s 

p— s 

* 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

r 

o 

o 

o 

H 

(VJ 

+ 

+ 

«n 

CO 

to 

to 

to 

N 

CM 

H 

o 

.1 

1  > 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

O 

<M  <M 

O  O 


ro  (\l  (M  (VJ 

o  o  o  o 

O  U  U  O 


2,  (MN(M(M(M(MNN<M(M 

^JOOOOOOOOOO 

zoouuouoooo 


I  tt  ti  ii  n  n  u  ti  it  it  n  it  it  it  n  it  it  n 


o 

o 

© 

p— «* 

pM» 

o 

o 

p«*s 

o 

p-s. 

O 

o 

o 

w 

CM 

CM 

o 

o 

o 

CM 

o 

"(M 

'nj 

0 

0 

o 

© 

CM 

M 

to 

CM 

CM 

0 

4> 

s 

s 

O 

2 

2 

u 

o 

u 

CM 

o 

2 

2 

2 

2 

2 

U 

2 

2 

2 

O 

4- 

4- 

4- 

4- 

4- 

o 

6 

4- 

4- 

4- 

4- 

4- 

4- 

4- 

4- 

4- 

4- 

pH*, 

pp*. 

p«» 

p-v 

o 

y 

pa* 

p_* 

.pa* 

PS 

p-N, 

p-S 

P-S. 

P-S 

H 

H 

o 

o 

H 

H 

H 

H 

H 

O 

o 

o 

o 

o 

O 

o 

o 

(\1 

m 

4- 

4- 

o 

o 

o 

o 

o 

o 

CO 

«M 

^4 

co 

o 

o 

o 

O 

o 

— 

o 

o 

o 

o 

o 

pH 

o 

O 

o 

o 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

O 

O 

O 

U 

O 

U 

u 

o 

U 

2 

U 

U 

o 

u 

u 

O 

o 

V 

o 

u 

o 

H 

CM 

m 

'p 

in 

vO 

«H 

N 

m 

m 

vO 

00 

O' 

pH 

pH 

H 

pH 

pH 

pH 

H 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

A  -14 


F 


t 

!  f 

I  *» 

< 

j. 

i 

i  i 


/.* 


r 


i. 


r. 

s  I: 

I  0 

L 

» 

I 

L 

L 
L 

r  • 
i 
J 

r 


4> 

OJ 


M> 

00 


m 

00 


oo 

r- 


O' 


o 

00 


1-1  <\] 
00  00 


m 

00 


o 

00 


o  M* 
00  00 


r-  o 
oo  co 


© 

00 


cs 

V 


Vi  O 


u 

o 

o 

V 

5 


<M 

I 


H 

\ 

o 

00 

in 

*w* 

8* 

o> 


O  NO  (M 

co  cm  co 

I  I  I 

00  H  H 

•  •  • 

CM  s0  (M 


H 

Pi 

> 

o 

r- 


S4 

4) 

m  i-H  >-i 

mm  m 

i  i  ■ 

o  in  a" 

•  •  • 

mm  co 


H 

OJ 

\ 

o 

o 

o 


N 

Ctf 

\ 

o 

o 

o 


H 

« 

o 

00 

00 


5  “1  v 

S'  ^  S' 

4)  H  4) 

M*  •-<  m 

m  mm 

i  i  i 


H 
Pi 
\ 
o 
o 

2  o 


o 

00 


H 

QiS 

N. 

O 

T 

m 

i 

Sm* 

S' 

4) 


O 

00 


*  S' 

fy  m 


cm  m  co 

mm  m 

■  i  i 


m 

oo 


00  o 

ft  ft 

.-t  (M 


©  o 

ft  ft 

-H  TJ. 


oo 


o 

m 


oo 

m 


>o 

2 

M3 

2 

+ 

sD 

2 

+ 

O 

2  * 

*  +„ 

M3 

2 

+ 

M> 

2 

+ 

+ 

CM 

o 

(M 

(M 

o 

§  1 

0 

rH 

X 

X 

U  X 

S 

£ 

11 

n 

11 

tin  n 

tl 

+ 

o 

(M 

O 

O 

2 

m> 

2 

+ 

M3 

2 

+ 

M> 

2 

+ 

CM 

CO 

+ 

(M  0  O  UO  O  <m 

a  u  o  a  a  o 


u 

4) 

ptj 


M3 

M3 

2 

r'- 

2 

2 

2 

M3 

M3 

M3 

M3 

5 

+ 

+ 

4* 

(M 

M3 

M3 

+ 

M3 

M3 

2 

2 

M3 

2 

+ 

+ 

a 

+ 

CM 

a 

+ 

o 

z 

0 

z 

2 

+ 

2 

+ 

o 

O 

2 

+ 

+ 

0 

4* 

CM 

o 

2 

+ 

a 

o 

o 

o 

+ 

+ 

z 

0 

O 

a 

z 

a 

o 

+ 

+ 

+ 

o 

CM 

o 

CM 

o 

+ 

+ 

+ 

+ 

+ 

•+• 

+ 

a 

a 

a 

U 

a 

a 

z 

z 

o 

0 

o 

0 

0 

o 

H 

CM 

CO 

H 

(M 

m 

in 

M3 

00 

O' 

H 

H 

H 

H 

►“J 

h 

t-5 

*“3 

•-3 

t-1 

*n 

A-16 


J 


J14  O  +  O,  +  M6  —  O,  +  M,  1.0(-34)  exp(1010/RT) 


a:  0000000040 

oo  oo  a*  00  00  O'  oo 


O'roCT'O'co^iininfoo 

N  0s  O'  O'  0s  O'  CO  00 


H 

H 

H 

H 

H 

H 

H 

H 

-*-> 

01 

E 

E 

H 

E 

E 

E 

e 

0) 

N 

N 

\ 

N 

E 

\ 

\ 

\ 

\ 

•>4 

o 

o 

o 

o 

\ 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

•H 

44 

o 

oo 

sO 

m 

o 

O' 

o 

o 

m 

m 

m 

pH 

pH 

r~ 

H 

pH 

pH 

04 

4> 

co 

1 

m 

1 

1 

i 

1 

1 

1 

■ 

o 

i 

>»P 

•mP 

w 

«w># 

U 

¥ 

(D 

S' 

o 

S' 

a. 

X 

S' 

a 

X 

¥ 

S' 

¥ 

<u 

0) 

H 

t) 

cj 

<u 

V 

V 

4) 

4-> 

-—V 

, — * 

^  . 

, _ „ 

p-s. 

. _ _ 

#■*>* 

«>— ** 

H 

CO 

CO 

o 

o 

H 

O  pH  pH 

o 

cO 

pH 

pH 

pH 

CO 

P5 

00 

1 

H 

•H 

H 

H 

H 

»H  pH  rH 

1  a  • 

H 

pH 

pH 

H 

pH 

pH 

J. 

, _ . 

w- 

■  II 

1 

-V 

1 

1 

1 

1 

irT 

co 

nO 

o 

p—4 

H 

CO 

o  vo  ro 

?T 

o 

cT 

o 

CO* 

o 

• 

• 

• 

• 

• 

•  • 

• 

• 

• 

• 

• 

cO 

CO 

»H 

o 

CM 

<M 

00 

vO  04  tT 

pH 

00 

00 

pH 

04 

o 

E 

04 

o 

E 

+ 

04 

E 

E 

O 

(M 

E 

O 

+ 

E 

+ 

E 

+ 

+ 

+ 

+ 

_ 

o 

O 

o 

o 

O 

o 

o 

O 

E 

E 

E 

o 

O 

0 

U 

a 

U 

U 

o 

o 

O 

0  N  ^  O 

2  O  o  O  + 


m 

5  2 


E  E  <M  E  O 
O  O  E  O  2 


04  O 


<M  O 

<M  O  nl  • 

w  /\  T 


E  (M  E  <M 
0  2  0  2 


li  II  tl  H  IHI  il  II  ll  t|  ll  ll  II  II  !|  tl  II 


f 

1  j 

.  2 

E 

+ 

0 

+ 

HH 

0 

+ 

0 

E 

0 

04 

o 

04 

o 

E 

! 1 

L  o 

C4 

O 

04 

0 

04 

o 

04 

+ 

0 

+ 

0 

+ 

o 

t 

E 

E 

E 

E 

E 

E 

E 

+ 

f 

j 

L 

U 

U 

U 

U 

O 

u 

E 

( 

*H 

04 

m 

in 

v£> 

0- 

oo 

1  E 

E 

E 

E 

E 

E 

E 

E 

in 

~  O 

O  (M  >r  m 

E  S  2  2  O  g  0 

+  +  +  +  +  +  + 

N  (M  (M  (M  (M  N  (M 

o  o  o  o  o  o  o 


00  O'  r— l  r-l  *— «  r-l  *— <  r-*l  w—i 

EEEEEEEEEE 


t  z1'^  » 


A  - 17 


V* 


I 

[ 


r 


i 


1  /  ^ 


rvj 

°  g 

+ 


**  o 
o  z 


£ 


o 

£ 


CO 

o 


CO 

o 

+• 

(M 

CO  O 

O  2 
2  + 

+  CO 
CO  O 

O  2 


£ 


CO 


o 

o 

o 


X 


CO 

£ 


fM  O 

o  u 


r1  o  o 


CO  (M 

£  X 


rO 


°c  o 


£ 


CO 

o 

£ 

£ 


(M 

r  o 


o 

+ 


CO 


£  £ 


a 

o 


s  n  ti  n  n  u  n  n  ii  ti  n  u  ti  u  u  n  u  ti 


i> 

pi 


I 

<N1 

g 

CO 

0 

£ 

(M 

o 

ro 

O 

CO 

0 

£ 

CO 

0 

£ 

in 

0 

CM 

£ 

CO 

0 

o 

0 

O 

+ 

CO 

£ 

+ 

co 

0 

+ 

CO 

o 

+ 

£ 

+ 

O 

£ 

+ 

CO 

O 

£ 

+ 

O 

+ 

£ 

O 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

fj 

£ 

£ 

£ 

£ 

0 

0 

0 

0 

0 

0 

0 

0 

O 

O 

O 

O 

O 

f 

00 

O' 

o 

CO 

CO 

'j' 

in 

o 

00 

O 

o 

co 

CO 

NT 

i 

H 

H 

CO 

CO 

CO 

CM 

CO 

IN] 

co 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

£ 

A  - 18 


s  v 


to 

o 


o 

oo 


© 

oo 


O' 

f- 


CO 

o  o 

00 


o 

00 


H 


H 

ptj 


H  H  H 
P5  PJ  BJ 


o 

o 

o 

<M 


O 

O 

M1 

<M 

t 


o  o 
o  o 

<M  O' 
CO  't 

I  I 


o 

o 

o 

sO 

I 


a 

5 


8-  S-  S' 


0)  0) 


0) 


rsj 


CO 


fM 

H  H  H 

I  I  I 

o  r-  ST 

•  •  • 

»-«  CO  C7^ 


CM 

H 

I 

o 


co  co  ro 

H  H  H 

I  I  I 

CO  oT  o 


(NJ  H 


m 


o 

z 


fM 


o 

2 


fM 

fM  O  CM 
0  2  0 


fM  fM  fM 

o  o  o 


+  +  +  in  +  +  + 

CM  fM  fM  O  fM  CO  (M 
O  O  O  CM  o  o  o 

x  z  z  z  z  z  z 


II  II  II  II  II  II  II 


CO 


CO 

o 


o 

z 


CO 

o 

J?  z 


CO 

§ 


CO  0 
0  2 


X 

o 


§ 


g 


fM 

g 


00 

o 

z 


ro 

g 


CO 

§ 


in 


vo  oo  O'  ©  i-i 

co  fO  rO  CO  co  M"  M1 
X  * 


1 


A-19 
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